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Two-loop anomalous dimensions and one-loop renormalization scheme matching factors are calcu¬ 
lated for six-quark operators responsible for neutron-antineutron transitions. When combined with 
lattice QCD determinations of the matrix elements of these operators, our results can be used to 
reliably predict the neutron-antineutron vacuum transition time, T„fi, in terms of basic parameters 
of baryon-number violating beyond-the-Standard-Model theories. The operators are classified by 
their chiral transformation properties, and a basis in which there is no operator mixing due to strong 
interactions is identified. Operator projectors that are required for non-perturbative renormalization 
of the corresponding lattice QCD six-quark operator matrix elements are constructed. A complete 
calculation of Sm = l/r„s in a particular beyond-the-Standard-Model theory is presented as an 
example to demonstrate how operator renormalization and results from lattice QCD are combined 
with experimental bounds on Sm to constrain the scale of new baryon-number violating physics. At 
the present computationally accessible lattice QCD matching scale of ~ 2 GeV, the next-to-next- 
to-leading-order effects calculated in this work correct the leading-order plus next-to-leading-order 
Sm predictions of beyond-the-Standard-Model theories by < 26%. Next-to-next-to-next-to-leading- 
order effects provide additional unknown corrections to predictions of Sm that are estimated to be 
< 7%. 

PACS numbers: 11.30.Fs,12.38.Bx 


I. INTRODUCTION 

The universe contains many more baryons than antibaryons [T]. Unless this baryon asymmetry is attributed to 
fine-tuning of the initial conditions of the universe, the baryon asymmetry must have been generated dynamically 
during the early universe. Any mechanism describing this process of baryogenesis must include violation of baryon- 
number (B) conservation, violation of C and CP, and departure from thermal equilibrium [2]. The Standard Model 
includes B violation through non-perturbative electroweak processes that violate B L but preserve B — L [ 3111 . 
It also includes classical C and CP violation and departure from thermal equilibrium during the electroweak phase 
transition. However, the B and CP violating effects present in the Standard Model cannot reproduce the observed 
magnitude of the baryon asymmetry m\- As a result, Beyond-the-Standard-Model (BSM) physics is needed to 
explain baryogenesis. BSM baryon number violation could occur in many different ways. Theories that allow AB = 1 
transitions can allow B — L conserving proton decay, ^ which has been experimentally constrained to a high degree P- 
E]. Other classes of BSM theories do not allow proton decay, but do allow other baryon number violating processes. 
These models often instead include the AB = 2, B — L violating, neutron-antineutron transition [T2l-|39]. 

In vacuum, neutron-antineutron (nn) transitions would manifest themselves as oscillations between neutrons and 
antineutrons. The probability that a free neutron has transformed into an antineutron after time t is given by 
Pnn = sTc? (t/Tnri), where Tnn is the neutron-antineutron vacuum transition time. Experimental measurements of 
magnetically shielded cold neutron beams at the Institut Laue-Langevin (ILL) have established a limit of Tnn > 
2.7 years [40) . There are also experimental bounds on the decay rate of neutrons bound in nuclei from large volume 
underground detectors. Super-K has bounded the transition time roie for nn transitions in oxygen, r^ie > 1.89 x 
10^^ years [41]. Nuclear structure calculations can be used to relate this nuclear transition time to the vacuum 
transition time Tnn- This bound on the vacuum transition time is estimated to be a factor of four or five larger than 
the ILL bound, but the nuclear structure calculations introduce non-trivial systematic uncertainties.^ It is believed 
that improvements in neutron transport/optics and neutron moderation technologies since the 1994 ILL experiment 
would allow for new neutron beam experiments to improve the ILL bounds by an order of magnitude or more [48] . 
There has been a recent push from both theoretical and experimental communities in support of new, state-of-the-art 
nn experiments. |43H45j . 

In order to constrain BSM theories predicting nn transitions, experimental results must be compared to theoretical 


^ See Ref. [8] for a recent review on proton decay 

^ In particular Ref. m cites a derived bound of > 7.7 years. More recent structure calculations in Ref. m modify this bound to be 
'^nn > 10.9 years, as noted in Ref. m- 
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predictions for Tnn- Making reliable predictions for Tnn within a particular BSM theory is challenging. In particular, 
theoretical descriptions of the nn transition process must include strong interaction physics as well as BSM physics. 
These effects are important at very different scales. High-scale BSM physics gives rise to effectively local AB — 2 
interactions turning three quarks into three antiquarks. Comparatively low-scale strong interactions bind these quarks 
(antiquarks) into a neutron (antineutron). Theoretical descriptions of this high- and low-scale physics can be factorized 
by using a Standard Model effective field theory description of nn transitions. In this approach, the Hamiltonian 
governing nn transitions is described as a linear combination of operators built from Standard Model fields. 

The most relevant Standard Model effective field theory operators contributing to nn transitions are dimension- 
nine six-quark operators. A complete basis of these six-quark operators can be constructed without specializing to 
a particular BSM theory. This construction was begun in Refs. [461117], generalized and detailed in Refs. [48l l4^ . 
and completed in Ref. |50| where spin-color Fierz identities were used to remove redundant operators from the 
basis. Higher-order operators of potential interest have also been discussed |CT1|S2|. The effects of low-scale strong 
interaction physics on nn transitions are encoded in quantum chrodynamic (QCD) matrix elements of six-quark 
operators between initial neutron and final antineutron states. All high-scale physics and BSM model dependence 
is encoded in the particular numerical coefficients used to express the effective Hamiltonian for a given theory in a 
six-quark operator basis. These numerical coefficients can be calculated perturbatively in BSM matching calculations 
for particular theories of interest. 

Testable predictions for Tnn cannot be made without reliable calculations of six-quark operator QCD matrix el¬ 
ements. Equivalently, experimental bounds on Tnn cannot be used to constrain BSM theory without reliable QCD 
matrix element calculations. These six-quark matrix elements have been estimated in the MIT bag model [151133]) but 
model estimates introduce uncontrolled uncertainties into the relation between BSM parameters and experimental 
observables |53| . The only available method to determine hadronic matrix elements with controlled uncertainties is 
lattice QCD. Preliminary lattice QCD calculations of nn matrix elements are underway [53]. Once completed, lattice 
QCD nn matrix elements can be non-perturbatively renormalized and then combined with BSM matching calculations 
performed with renormalized perturbation theory. 

The need for perturbative nn operator renormalization arises because lattice QCD matrix elements can only be 
renormalized at scales smaller than the UV cutoff of the lattice, (typical calculations today use lattice matching scales 
of po — 2 GeV [55]) but renormalization scales that are currently accessible in lattice QCD simulations cannot be 
(directly) used for perturbative BSM matching calculations. These BSM matching calculations receive logarithmic 
corrections that become large enough to invalidate perturbation theory unless the renormalization scale chosen is 
comparable to high scales where BSM physics becomes important. For typical BSM theories, these scales are in the 
range Absm = 10^ — 10^® GeV. To address this issue, renormalization group (RG) techniques can be used to sum 
these large logs and reliably relate matrix elements calculated with different renormalization scales. This RG evo¬ 
lution (“running”) and typical BSM matching calculations are both simplest in mass-independent, renormalization 
schemes such as modified minimal subtraction (NDR-MS).® The MS renormalization scheme can only be applied 
directly to dimensionally regularized matrix elements, and in particular cannot be applied directly to lattice regu¬ 
larized matrix elements. Instead, the Regularization-Independent-Momentum (RI-MOM) scheme can be introduced 
as an intermediate renormalization scheme |56| . As long as the lattice matching scale po used for non-perturbative 
renormalization is larger than hadronic scales where QGD becomes non-perturbative, it is possible to relate RI-MOM 
and MS renormalized matrix elements perturbatively (“matching”). The perturbative calculation of RG running and 
matching factors therefore allows non-perturbatively renormalized lattice QCD matrix elements to be combined with 
perturbative BSM matching calculations to provide testable predictions for Tnn in BSM theories of interest. 

The largest corrections to t„h arising from RG evolution are encoded in perturbative one-loop-running factors. 
These have been correctly calculated for nn operators in Ref. |50j . One-loop running provides an overall multiplicative 
correction to non-perturbatively renormalized matrix elements, see Eq. Further RG corrections to this result can 
be organized as a power series in as(po)- In order to verify that this perturbative expansion is well-controlled at a 
given po, it is necessary to determine the first term in this as(po) power series. This term is parametrically O(as(po)), 
and includes one-loop-matching effects. When running to high scales /i where as(/r) <i; as(po), two-loop-running 
effects also contribute at O(as(po)) and must be included as well, see Eq. (H). This work provides the first calculation 
of the one-loop-matching and two-loop-running factors needed to reliably estimate the convergence of RG relations 
between nn matrix elements at low scales po accessible to lattice QCD simulations and high scales p accessible to 
perturbative BSM matching calculations. 


® Naive dimensional regularization (NDR) prescribes that 75 anticommutes with iri D dimensions. Since closed fermion loops do 
not appear in nn calculations, no complications arise from using the NDR prescription. In the remainder of this paper we abbreviate 
NDR-MS as MS for brevity. 
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The remainder of this paper begins with a summary of our final results in Sec. [TTj Results are presented for the 
fixed-flavor basis commonly used in the literature and for a new chiral basis that is diagonal under RG evolution. 
The construction of this chiral basis is presented in Sec. Ell The RI-MOM renormalization scheme and associated 
operator projectors needed for perturbative and non-perturbative nn operator renormalization are defined in Sec. |W 
Calculation of one-loop-matching factors relating RI-MOM and MS renormalized operators is discussed in Sec.^ 
Calculation of two-loop-running factors is discussed in Sec. |VlJ Both Sec. [V|VI| discuss the careful treatment of 
evanescent operators vanishing in = 4 that is necessary for a correct calculation of RG effects. To demonstrate 
the phenomenological application of our results, a complete calculation of and resulting experimental constraints 


are discussed for a simplified BSM model in Sec. VII Physical results and implications are summarized in Sec. VIII 


Analogous one-loop-matching and two-loop-running calculations have been performed for four-quark weak matrix 
elements [57H55] and proton decay [5SH75] . These provide useful techniques as well as cross-checks for intermediate 
results. We avoid discussion of established techniques for multi-loop diagram evaluation in the main text, but for 
readers unfamiliar with two-loop diagram evaluation we present a pedagogical discussion in Appendices 0-0 Our 
explicit evanescent operator basis (technically required for a full definition of MS operator renormalization) is presented 
in Appendix]^ and some intermediate results are shown in Appendix [P] 


II. SUMMARY OF RESULTS 


The neutron-antineutron vacuum transition time predicted by a particular BSM theory can be calculated from 
matrix elements of the Hamiltonian density 




( 1 ) 


where the Qi (/r) form a complete basis of dimension-nine local six-quark operators with non-vanishing matrix elements 
{n\Qi{^)\n) between initial neutron and final antineutron states, the Ci{^) are Wilson coefficients, and /r is a 
renormalization scale. The Wilson coefficients are renormalization scheme and scale dependent and will differ between 
BSM theories. They can be calculated by matching tree- or one-loop- level nn amplitudes between the full BSM theory 
and an effective theory containing only Standard Model degrees of freedom. Hadronic matrix elements of are 

independent of the BSM theory used to calculate the C'/(^) but renormalization scheme and scale dependent. 

Lattice QCD first determines matrix elements of bare, lattice regularized operators. By subsequent lattice QCD 
calculations, these bare matrix elements can be non-perturbatively renormalized in the RI-MOM scheme described 
in Sec. IV at a lattice matching scale po- Provided as{po) ‘C 1, dimensionally regularized perturbation theory can 
be used to relate RI-MOM renormalized matrix elements to MS renormalized matrix elements. Introduction of RI- 
MOM as an intermediate renormalization scheme is necessary because the MS scheme can only be directly applied to 
dimensionally regularized (and not, for instance, lattice regularized) matrix elements. Setting the MS renormalization 
scale p = po removes large logarithms from the RI-MOM matching calculation. Perturbative calculations of 
in a particular BSM theory typically introduce additional logarithmic corrections hi{p/K bsm)- Since lattice QCD 
computational limits demand po ^bsm, Wilson coefficients calculated at p, = ^bsm must be RG evolved to p = po 
and then combined with MS renormalized matrix elements to include all large logs in BSM theory predictions of r„ 7 j. 

The renormalization scale dependence of the Wilson coefficients is encoded in the MS anomalous dimension matrix 
77 j, defined in Sec. lYll In Sec. m we use chiral flavor symmetry to construct an operator basis where the anomalous 
dimension matrix is diagonal. The RG equations relating Wilson coefficients at different renormalization scales can 
be solved perturbatively in this diagonal chiral basis. Including one-loop-matching and two-loop-running effects, the 
relation between the desired Hamiltonian the BSM matching coefficients at arbitrary scale p, and the 

non-perturbatively renormalized operators Q^^{po) used in lattice QCD simulations is 
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Chiral Basis 

Flavor Basis 

( 0 ) 

Ti 


^( 0 ) 

Qi 

^RRRy ^LLL 

4 

335/3- 34N//9 

101/30+ 8/15 In 2 

Q2 

0\rr, Orrr, Orrr, Orrr 

-4 

91/3-26/V//9 

-31/6+ 88/15 In 2 

Qs 

0\lr^ ^RRL 

0 

64- 10iV//3 

-9/10 + 16/5 ln2 

Qi 

( 4/5 Orrr + 1/5 Orrr) , 
( 4/5 Orrr + 1/5 Orrr) 

24 

229 - 46Nf /3 

177/10- 64/5 In 2 

Qs 

Orrr, Orrr, Orrr, 

Orrr, Orrr, Orrr, 

(2/3 Orrr + 1/3 Orrr), 

(2/3 Orrr + 1/3 Orrr), 

(2/3 Orrr + 1/3 Orrr), 

(2/3 Orrr -|- 1/3 Orrr) 

12 

238 - UNf 

49/10 - 24/5 In 2 

Qi 

(1/3 Orrr — 1/3 Orrr), 

(1/3 Orrr — 1/3 Orrr) 

4 

797/3 - 118Nf/9 

-109/30 +8/15 ln2 

Q3 

(1/3 Orrr — 1/3 Orrr), 

(1/3 0%RR — 1/3 Orrr) 

0 

218-38N//3 

-79/10+ 16/5 In 2 


TABLE I: Summary of results. The left-most column lists the chiral basis operators Qi with independent NNLO operator 
renormalization factors. The second column lists the corresponding Hxed-fiavor basis operators used in Ref. | 46H5U| that 
renormalize identically to Qi, see Sec. Ill Each Qi is equal to (—4) times the first fixed-flavor basis operator listed in the 
corresponding row of the second column. The other fixed-flavor basis operators listed may not be directly proportional to Qi 
but share the same one-loop MS anomalous dimension 7 }°^ (third column), two-loop MS anomalous dimension 7 }^^ (fourth 
column), and one-loop Landau gauge RI-MOM matching factor (fifth column) appearing in Eq. §. and depend 
on the evanescent operator basis used to extend D = 4 Fierz relations to D-dimensional operator relations in dimensional 
regularization. Our evanescent operator basis is presented in Appendix One-loop BSM matching calculations must use 
the same evanescent operator basis for consistency. Tree-level BSM matching calculations are unaffected, see Sec. m for an 
example matching calculation. 


where is a one-loop-matching factor defined in Sec. 7 !*^^ and 7 ^^^ are one-loop- and two-loop-running factors 
defined in Sec. lYll and /3o and /3i are well-known perturbative coefficients of the QCD /3-function presented for 
reference in Eq. (46). Only /3o, and 7 ^^^ depend on the number of active quark flavors, Nf. Matching between 
theories with different Nf at quark thresholds is included in the same manner as in RG evolution of weak matrix 
elements without penguin contributions since no penguin diagrams exist for nn operators. 

Ignoring QCD effects on RG evolution gives the leading-order (LO) result Ui{fi,po) = 1. Next-to-leading-order 
(NLO) QCD effects give a multiplicative correction to Ui{p,po) whose size is determined by the one-loop-running factor 
7 !*^^ correctly calculated in Ref. [50] . Higher-order corrections due to matching and running provide additive corrections 
that can be perturbatively expanded in powers of 08 (^ 0 ) and as{p). For high scales p where as{p) ^ asipo), 
Eq. (§ shows that one-loop-matching and two-loop-running effects receive similar 0 (q!8(po)) suppression. Both one- 
loop-matching and two-loop-running effects must therefore be included in a next-to-next-to-leading-order (NNLO) 
calculation of Ui{p,po). Next-to-next-to-next-to-leading-order (N^LO) corrections not included in Eq. ([^ arise from 
two-loop-matching and three-loop-running effects that are both O{as{po)^) suppressed. 

The NNLO operator renormalization factors and 7 }^^ are calculated for the first time here and summarized in 
TablelU The relative size of NNLO to NLO corrections to Ui{p,po) depends on p and differs between operators. Taking 
Po = 2 GeV and using the four-loop parametrization of as{p) in Ref. [75], NNLO corrections to NL O-fLO results 
for 6m = IjTnn are < 26% for all p > Po and may be significantly smaller in some BSM theories. Sec. VII presents 
an example calculation of the nfi vacuum transition rate for one of the simplified models of Ref. [25] . In this model 
the relative size of NNLO to NLO corrections to Sm is 14%. Estimating that unknown N^LO O(as(po)^) corrections 
are comparable to the square of NNLO O{as{po)) corrections allows systematic uncer taint y due to unknown N^LO 
corrections to be quantified as < 7% generically and 2% in the model discussed in Sec. IVII] 
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III. CHIRAL OPERATOR BASIS 


The operators relevant for nn transitions are Lorentz, color, and electromagnetic singlet six-quark operators of 
dimension nine. Since hadronic matrix elements must be calculated in lattice QCD simulations that only maintain 
approximate chiral symmetry at best, operators that are not singlets of the full electroweak gauge group should be 
considered. Even so, operator renormalization is most simply performed in the limit of massless up and down quarks. 
Classifying operators according to the SU{2)l x SU{2)fi chiral symmetry of QCD in this limit proves quite useful.^ 
In this section we construct a basis of irreducible chiral tensor operators that do not mix under perturbative QCD 
interactions. Fierz relations and symmetries of the color, spin, and flavor tensors used throughout this section are 
detailed in Appendix [X} Our notational conventions are as follows: we use as fundamental color indices, 

p,... as 4- vector Lorentz indices, a, /3,j,... as Lorentz spinor indices, a,b,c,. .. as flavor spinor indices, 
as operator basis labels, and y’s as chirality labels L, R. We will use A,B,C,... to denote adjoint indices in both color 
and flavor. su(3)c color generators will be denoted by and normalized to Tr(t^t'®) = while su{2)l and su(2)i^ 
flavor generators will be denoted by and normalized as Pauli matrices = 25^^ with ± ir^). 

We use Euclidean 1-) metric signature and will not distinguish between raised and lowered indices. Final results 

are valid in Minkowski signature; intermediate steps are not. Summation convention applies to all indices but not to 
operator basis /, J, AT,... and chirality y labels. 

Two quarks can be combined into a spin-singlet diquark by contraction with the antisymmetric charge conjugation 
matrix C and projected onto definite chirality by including Pl,r = ^(1 T Ts)- In D = 4, spin Fierz identities can 
be used to express any product of vector diquarks containing 7 ^ or tensor diquarks containing as a 

product of scalar diquarks. Denoting flavor doublet quark fields by = (uf,df), only operators containing three 
products of scalar diquarks need to be considered. 

Flavor Fierz identities allow us to only consider operators where each diquark is either a flavor singlet contracted 
with the antisymmetric tensor or a flavor vector contracted with the symmetric tensor 


( 3 ) 

where we have suppressed free color indices. Irreducible su(2)^-spin-two and su(2)^-spin-three chiral tensor operators 
can then be dehned as 


V 


AB 

X 


A, /K Q 


C 
X ’ 


j^ABC _ 

— X X X 


SAB^{C^D^D} 

A A- A- AAA A A A 


( 4 ) 


Since operators contributing to nn transitions must lower the third SU{2)v isospin component isospin by one 
unit,® at least one diquark must be contracted with to form a dfd^ diquark. The other two diquarks 

must combine to have no net effect on I 3 . Taking this dfd^ combination to be our third diquark for convenience 
and enforcing antisymmetry under quark exchange, the only available tensors for constructing color singlet six-quark 
operators are 


rpSSS 

^{ij}{kl}{mn} 

rpAAS 
^ [ij][kl]{mn} 


— ^ikm^jln H” ^jkm^iln ^ilm^jkn “t“ ^ikn^jlnm 

— ^ijm^kln ^ijn^klmi 


( 5 ) 


where { } denotes index symmetrization and [ ] denotes index antisymmetrization. From here onward we suppress 
explicit quark indices and use the diquark notation 


{i^.CPrit%) = mCPRr^[^rXt^^^,. 


( 6 ) 


We further suppress color indices in diquark products, e.g. 


We thank Brian Tiburzi for very helpful insights on these chiral transformation properties. 

^ In particular, nn transitions only involve operators with negative parity, A/ = 1, and A/a = — 1. We only explicitly enforce the latter 
constraint A /3 = — 1 in order to simplify the perturbative calculations presented here. 
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Chiral Basis 

Fixed-Flavor Basis 

Chiral Tensor Structure 

Chiral Irrep 

Qi 

^RRR 

VrVrV+T^^^ 

(1l, 3h) 

Q 2 

^LRR 

VrVrV+T^^^ 

(1l, 3h) 

Qa 

0\lr 

VrVrV+T'^^^ 

(1l, 3h) 

Qi 

4/5 Orrr 4- 1/5 Orrr 

'j^33-\- rjnS S S 

{1e,7r) 

Qs 

(^RLL 


(5i,3H) 

Qe 

Orel 


(5i, 3 _r) 

Q 7 

2/3 Orer -b 1/3 Oeer 


(5e,3r) 

Qi 

1/3 Orrr — 1/3 Orrr 

VrVrV^T^^^ 

(1e,3r) 

Qa 

1/3 Oeer ~ 1/3 Orer 

VrVrV+T^^^ 

(1e,3r) 


TABLE II: The chiral basis operators Qi shown in the first column are equal to (—4) times the corresponding fixed-fiavor 
basis operator combinations shown in the second column. Each chiral basis operator is equal to a color contraction of the 
tensor operators ' shown in the third column. The corresponding chiral irrep of each operator is shown in the last column. 
Qir ■ ■ ) Qi a-nd their parity conjugates (L •f-J- i?) form a complete basis for nn transition operators in D = 4. Since they are 
components of the same chiral tensor operator , Qe and Q 7 renormalize identically to Q 5 and are redundant for our 

purposes. Qi and Qa are equal to Qi and Q 3 in D — 4, but they renormalize independently in MS at NNLO. 


Using these building blocks and neglecting operators that have A /3 — 1 or vanish by quark anticommutivity, we 

find that at NLO there are five chiral tensor operators with independent renormalization properties, 


Qi = {ipC , (7a) 

Q 2 = {'ipCPLiT‘^^l)){'il)CPRiT‘^'4)){'4>CPRiT‘^T+'il))T^^’^, (7b) 

Qz = {iiC PLiT'^ij^){'il}C PLiT‘^'ip){'ipC Prit'^t^ , (7c) 

Qi = {'tpCPRiT‘^T^ip){ipCPRiT‘^T^ilj){ipCPRiT'^T~^ip)T^^''^ (7d) 

- ^{'ipCPRiT'^T^'lp){'lpCPRiT^T^i’){'ipCPRiT‘^T~''^p)T^^^, 

Q 5 = {'tpCPRiT'^T~'ijj){ipCPLiT'^T^'4’){'ipCPLiT‘^T~^tp)T^^^. (7e) 


Symmetries of and under diquark exchange ensure that all products of flavor vector diquarks are totally 

symmetric. Qi includes a flavor trace subtraction. This ensures that all operators are irreducible chiral tensor 
operators. 

There are two additional operators that cannot be expressed as linear combinations of Qi, • ■ •, Qs, 

gg = {i^CPRiT‘^T^ij){i^CPLiT‘^T^ij){i;CPLiT^T+'tp)T^^^, ( 8 a) 

Qy = {lpCPLiT‘^T^'ip){'ipCPLiT‘^T^tp){tpCPRiT‘^T'^'tp)T^^^ ( 8 b) 

- \{i^CPLiT^T^-tp){ijCPLiT^T^i;){'iljCPRiT^T+i;)T^^^. 

O 

These two operators and Qg are different components of the same chiral tensor operator . This implies that 

Qs) Qsj Q 7 have identical anomalous dimensions and matching factors in renormalization schemes respecting 
chiral symmetry. In /) = 4, Qi, • • • , Q 7 and their seven parity conjugates found by taking L ^ R everywhere and 
including a relative minus sign form a complete basis of dimension-nine operators contributing to nn transitions. 

We also consider two more operators Qi and Q 3 that in 14 = 4 are equal to Qi and Q^ respectively by Fierz 
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identities, 

Qi = (9a) 

Qs = ^{'ipCPLiT‘^T^'ip){ipCPLiT^T^4’){i’CPR'iT'^T~''tp)T^^^. (9b) 

The Fierz relations Qi = Qi and Q 2 = Q 3 are broken in dimensional regularization, and Qi and Qa ^re independent 

operators in D dimensions. In principle, we could choose our physical operator basis to be Qi, ■ ■ ■ ,Q 5 and include 
Qi — Qi and Qs — Q 3 as additional evanescent operators vanishing in Z? = 4 but present in D dimensions. In practice, 
it is much easier to directly determine matrix elements of Qi and Qa and explicitly include them in the physical 
operator basis. For the purposes of NNLO operator renormalization we take our chiral basis operators Qj to include 

The basis commonly used in the literature involves fixed-flavor quark fields |46H5()j . 

= iuCP^,u){dCP^,d)idCP^,d)T^^^, 

^Lx.xa = {uCP^,d)iuCP^J){dCP^,d)T^^^, (10) 

^Lxaxa = iuCP^,d){uCP^,d)idCP^,d)T^^^. 

These fixed-flavor basis operators satisfy the relations = ^xRL ^lrx ~ ^RLx' ^ = 4, they also satisfy 

the Fierz identities , — 0\ > = These relations reduce the number of linearly independent operators 

to 14. It is straightforward to evaluate the flavor contractions oi ip = (u, d) in the Qi and verify that Qij ■ • ■ i Qt 
their parity conjugates form 14 linearly independent combinations of fixed-flavor basis operators. One can similarly 
verify that the Fierz relations Qi = Q\ and Q 3 = Q 3 are equivalent to the hxed-flavor basis Fierz relation above. The 
precise relations between the chiral basis and fixed-flavor basis operators and their explicit chiral tensor structures 
are shown in Table IE 


IV. RENORMALIZATION SCHEMES 


The commonly used MS renormalization scheme simplifies RG evolution, preserves important symmetries of chiral 
gauge theories such as the Standard Model, and is technically simple to implement in perturbative calculations 
performed with dimensional regularization. MS is limited, however, in that its defining renormalization condition 
can only be applied to regularized matrix elements calculated with dimensional regularization and not other with 
regularization schemes such as lattice. The RI-MOM renormalization scheme [S5], while not as technically simple to 
apply to dimensionally regularized matrix elements, has the advantage of a regularization independent renormalization 
condition. In this section we construct a RI-MOM operator renormalization condition for the Qi that can be applied 
both non-perturbatively to lattice QCD matrix elements and perturbatively to dimensionally regularized matrix 
elements. This is an essential intermediate step in connecting lattice regularized and MS renormalized nn matrix 
elements. 

In this section we explicitly display the spacetime and renormalization scale dependence of the quark fields ipf^{x, fP) 
and six-quark operators Qi{x,^). These renormalized quark fields and six-quark operators should be distinguished 
from their bare (regularized) counterparts, defined for quark fields of flavor q = u^dhy 


= Zg ^/2(/r)[q°]“(a;), Qiix,fi) Zij{n)Q°j{x), 

J 


( 11 ) 


where the wavefunction renormalization factor Zq{fj,) and operator renormalization matrix are formally de¬ 

fined by the renormalization conditions of a particular renormalization scheme. We denote perturbative expansion 
coefficients for either renormalization factor by 


Z{^) — I-\- 


47r 




sid) 


47r 


Z(2) -p 0(a3). 


The RI-MOM scheme wavefunction renormalization factor Z^^(/i) is defined by 


1 =—Tr ( 
48 ' ’ 


dSq{po,n=po) 


-1 


dp^ 


= ^Zf\po)Tr [ 7 /^^'^°^^°^ 


48 


dpi^ 


( 12 ) 


( 13 ) 






where Tr denotes a trace over color and spin (quark flavor q is held fixed), po is the lattice matching scale, and 
Sqipo) and Sq{po,p) are bare and renormalized quark propagators respectively. The normalization is chosen such 
that = 1 + 0{as). Eq. (13) assumes that a gauge fixing condition has been imposed so that the quark 

propagator is non-vanishing. In this work we consider a general gauge where the tree-level gluon propagator is 
5^^ {g^Lv/p^ ~ (1 ~ ClPfj.Pi^lP'^)- One-loop matching is performed in Landau gauge, ^ = 0. Two-loop-running results 
are gauge invariant, and for simplicity the two-loop calculation is performed in Feynman gauge, ^ = 1. 

A one-loop calculation of the quark self-energy in H = 4 — 2£ dimensions shows that the counterterm needed to 
renormalize the bare propagator according to the renormalization condition Eq. (131 is |56j 


« 3 



(14) 


where 1/e = 1/e —7£)-|-ln47r. The MS wavefunction renormalization factor is defined by the condition that Zq remove 
precisely the poles in 1/e from the quark propagator. At one-loop 


« 3 \e) 


(15) 


At two-loop order the quark propagator includes diagrams with divergent one-loop subdiagrams. These diagrams 
include non-local divergences proportional to ln(^^/p^)/e. Renormalizability guarantees that these non-local two- 
loop divergences cancel after including counterterm diagrams in which divergent subdiagrams are replaced by their 
one-loop counterterms The remaining local divergences are removed by a two-loop counterterm that in Feynman 

gauge is given by m 


vMS.(2) _ 44 -47 -I- 2Nf 

« 9e2 ^ 3e 


(16) 


where Nf is the number of active quark flavors. 

A regularization independent definition of Z|j(^) can be given in terms of a renormalization condition applied to 
vertex functions for each Qj. These vertex functions can be constructed, perturbatively or non-perturbatively, by Wick 
contracting Qj with interpolating operators for initial neutron and final antineutron states. A vertex function with 
Qj inserted at zero-momentum can be constructed by including three external antiquark fields carrying momentum 
p and three external antiquark fields momentum —p. These antiquark fields act as interpolating operators capable 
of creating initial neutron and final antineutron states. In order to simplify the non-perturbative construction of 
this vertex function in lattice QCD calculations, it is convenient to work with interpolating operators built from 
fixed-flavor quark fields. The quark fields must be assigned momenta symmetrically in order for the RI-MOM scheme 
defined using the vertex function to preserve chiral symmetry. A suitable definition is given by m 


= I (Qi{0)<{p)Ujip)dl{p)d^i{-p)(C{-p)di{-p)\ 

^ ' ai 

+ ^ {QiiO)uf{p)u^{-p)^kip)dUp)^m(.-p)dii-p)) 
+ ^ (Qi(0)u^(-p)u^(-p)dk(p)di(p)d^(p)di(-p)j 


(17) 


where antisymmetrization of all antiquark fields of the same flavor is implied, renormalization scale dependence is 
suppressed, and the subscript amp refers to the prescription of amputating external legs with the replacement 

g“(p,p) ^ 9“'(P,M)[5'7^(p,/i)]“'“. (18) 

When the Wick contractions for A/ are enumerated, the 20 distinct momentum routings available for Feynman 
diagrams with six indistinguishable external legs carrying momentum {p, p, p, —p, —p, —p} each appear with equal 
weight. Perturbative contributions to A/ are defined by 

Ai{p,p) =Af^ + a\^\p,p) + a\^\p,p). (19) 


® For a comprehensive review of renormalization theory with further references to the original literature, see Ref. Ea 
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It should also be noted that A/ is defined with an “exceptional momentum configuration” where the momenta of some 
subsets of external fields add to zero. However, vertex functions for purely baryonic operators like Qj are not expected 
to suffer from the non-perturbative chiral symmetry breaking artifacts that affect mesonic operators in exceptional 
momentum configurations. In particular, infrared divergences in the chiral limit arising from pseudo-Goldstone poles 
can lead to enhanced non-perturbative chiral symmetry violating mesonic operator mixing in lattice QCD simulations 
with chiral fermions, see for example |62j . Enhancements arise from diagrams in which an external quark and antiquark 
can be combined in a subdiagram with zero external momentum. For the purely baryonic operators considered here 
as well as for proton decay operators, there are no pseudo-Goldstone pole enhancements from diagrams in which two 
external quarks with positive baryon number are combined in a zero-momentum subdiagram. 

The RI-MOM scheme is defined by a renormalization condition on A/ (p, /i), 


SiJ =Tr ['PiAj{po,fi = po)] 

^J^ZfkiPo) (^“(po))'Tr [VjAlipo)] , 


( 20 ) 


K 


where A^(po) is a bare vertex function built from bare operators and amputated with bare propagators, Tr [ViAj] = 
Kfi^mn ’ Operator projectors Vj are defined by 


Tr 


ViA^P 


= Sij, 


( 21 ) 


Each external quark in Eq. (20) carries momentum ±po and the renormalization scale is identified with this lattice 


matching scale p = pg. As discussed above, po must be chosen to be much larger than hadronic scales to allow for 
perturbative matching but much smaller than the inverse lattice spacing used for non-perturbative renormalization 
to control discretization errors. For many quantities, these constraints are satisfied at po — 2 GeV. Gomparison of 
the size of 0{as{po)) NNLO corrections to the NLO result in Eq. Q should provide an estimate of perturbative 
convergence with a chosen pg. 


A set of projectors satisfying Eq. (21) for the chiral basis operators is given by 


(T> 

V' ^Hjklmn 

(T> \a0'rSri(^ 
V' ‘^/ijklmn 

(T> 

V ^Jijklmn 
(p> \al3'rSri(^ 

(T> \al3'rSri(^ 
\' ^ )ijklmn 

(T) 

V ^ ) ijklmn 


1 


92160 

1 

'18432 
1 

'36864 

1 

'221184 

1 

'221184 
1 


rpSSS 

^{ij}{kl}{mn} 


2Tr 


AAS 

ij][kl]{mn} 


icpRr\cpRr^{cpRr^ 


rpSSS 

^{ij}{kl}{mn} 


iVj) 


OiQ')5r)(^ 

ijklmn 


55296 

1 

'73728 


{CPRrP{CPRy\CPRT^ 

{CPLr\CPRy^{CPRr< + 2T[(^f;]{_j(CPL)“'(CPfl)'^«(CPfl 
-T[^}%inmn}iCPLr^iCPLr\cPR,r^ + 2r[(^f,]{_j(CPi)“'(CPi)'^^(GP^)’'^) , 
{^mkinmn}icpRr^icpRr\cpRr^ + sT^fg,^^^„^icpRr\cpRr^icpRri) , 
{Ti%%inrnn}iCPRr^iCPLrHCPLr^), 

{T{%%inran}iCPRr\CPLr^iCPLr^ + 6P[(^f„ (CPfi)“'(CP^ )'^^ (CP^)’''") , 


{T[%%i}{mn}iCPLr^iCPLr\CPRr^ 


2Pr 


AAS 

ij] [kl\{mn} 


{cPLr\cPLT^{cPRy^ 


(22a) 

(22b) 

(22c) 

(22d) 

(22e) 

(22f) 

(22g) 


We explicitly include projectors for Qg and Q^ since they must be analyzed separately in lattice QGD calculations 
without exact chiral symmetry. The seven parity conjugates of Qi,..., Qy should be analyzed separately in lattice 
QCD calculations; projectors for these operators are found by taking L ^ R everywhere and including a relative minus 
sign. Pi and Pg are suitable projectors for Qi and Qg since they are equal to Qi and Qg at tree level. Projectors 
and for fixed-flavor basis operators differ from those of Eq. (22) by an overall normalization factor of (—4) are are 
described in more detail in Ref. [50] . 

ZYf can be defined through a renormalization condition for dimensionally regularized vertex functions: at each 
order of renormalized perturbation theory, add counterterms that remove precisely the 1/e poles proportional to 
Aj from A/. A more precise definition of both the RI-MOM and MS renormalization conditions for dimensionally 
regularized amplitudes requires a careful treatment of evanescent operators. This is postponed to Sec. |VB| 
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V. ONE-LOOP MATCHING 


RI-MOM and MS renormalized operators with renormalization scale fj, = po are related by Eq. (|ll|) 

i;: 


QfHpo) = 

J,K 


Q^^{Po) = ^rijQ^^{po). 


(23) 


The matching factor rjj relates renormalized operators and is therefore a finite quantity, rjj can be consistently 
calculated perturbatively in terms of and Z^^ as long as both contain the same UV divergences and in particular 
are calculated with the same regularization. This allows us to perturbatively express rjj as 


rij{as) = 1 + 


= 1 + 


^ij 


ij 


asiPo) /^Ri,(i) _ 2 -ms,(i) 
47r 

asiPo )^ 

Att 


0(«s) 




(24) 


(0) 


Since the chiral basis operators do not mix under renormalization, Zjj and rjj are diagonal and we further define Tj 
to be the diagonal elements r^jj = (no summation on I). Differences between definitions of the renormalized 

coupling constant as{po) in different sche mes are formally 0{a^) and can therefore be neglected in Eq. (24|. When 


calculating numerical results in Sec. |VII[ we use an MS coupling constant definition for both two-loop running and 
one-loop matching. This defines the one-loop-matching factor appearing in Eq. (§ in terms of Zj^J. The 
remainder of this section describes the diagrammatic evaluation of zjy* from one-loop corrections to . 


A. Diagram Evaluation 

Eeynman diagrams representing corrections to Qi involve six quark lines carrying baryon number into a local 
vertex where the quark lines are joined to form three spin-singlet diquarks. To simplify the structure of these 
diagrams it is convenient to introduce charge conjugate quarks ('0‘^)fa = Expressing all diquarks in Qj 

_ Q 

as {ipCPxi’) = (V' Px'^) removes the need to introduce spin-transposed propagators and explicit factors of C at 
the six-quark vertex. With this approach, one quark line contained in each spin-singlet diquark is replaced with a 
conjugate-antiquark line with it’s fermion charge arrow pointing out of the six-quark vertex, as shown in Fig.[l] These 
obey standard Feynman rules for quark lines carrying fermion charge away from the vertex, except that conjugate 
quark-gluon vertices receive an extra minus sign and transposition of because 

Matching between RI-MOM and MS is performed in the limit of massless quarks where SU{2)l x SU{2)ji chiral 
symmetry guarantees that loop-level operator corrections will contain diquarks with the same chiral structure as 
the tree-level operator. When calculating diagrams with no closed fermion loops in NDR, explicit factors of 
can therefore always be moved from the operator vertex to one end of the quark line representing each spin-singlet 
diquark. In addition, the tree-level flavor structure of each operator is preserved diagram-by-diagram because of the 
flavor-blind nature of QCD. With these considerations, the only non-trivial tensor structure that needs to be inserted 
at each six-point operator vertex is (1 0 1 0 for Qi, Q 2 , and Q 3 and (1 (g) I 0 for Q 4 , Qs, Qi, and 

Q 3 . Diagrams with six-quark vertex factors of (1 (g) 1 (g) and (1 (g) 1 (g represent amplitudes denoted by 

and A4'® respectively. These amplitudes provide loop-level operator corrections to Qj once factors of CP^, flavor 
tensors, and contractions with quark fields are included. As a concrete example, diagrammatic operator corrections 
to Qi are found from the amplitude by making the replacement 

(Ti g Ta g TgjT ^ {-^pCPRTliT^P:){PJCPR^ 2 iT^^p){i^CPR^ 3 iT^T+PJ)T^^^. (25) 

To determine finite 0(e°) contributions to aP , all contributing diagrams should be calculated with all distinct mo¬ 
mentum routings assigning incoming momenta {p,p,p, —p, —p, —p} to the external legs. aP is then found by adding 
external quark fields to build the correlation functions of Eq. 0. performing Wick contractions, and amputating 
external legs. The contribution of each Wick contraction is represented by the sum of all contributing amputated 
diagrams with a particular momentum routing. 1/e pole terms are momentum independent, and can be determined 
from any momentum routing free of infrared divergences. 

The topologically distinct classes of one- and two-loop Feynman diagrams contributing to aP and aP are shown 
in Fig. Calculating rP requires evaluating the 15 one-loop diagrams in classes d = 1-3. For each of the one-gluon- 
exchange diagrams in d = 1-3, the two distinct momentum routings correspond to gluon exchange between quarks 
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LO 




NLO 


NNLO 





FIG. 1: The tree-level operator diagram, 15 one-loop diagrams, and 320 two-loop diagrams evaluated in this work. Operator 
insertions are represented by a six-point vertex joining three quarks and three conjugate antiquarks. The operator insertions 
are local; the separate solid lines represent propagators for quarks contracted into separate spin-singlet diquarks as indicated by 
the fermion charge arrows. The external quarks are assigned momenta {p,p,p, —p, —p, —p}. Diagrams are organized into classes 
that share the same loop integrals and Dirac structures. The number of diagrams in each class is shown. All two-loop diagrams 
with divergent subdiagrams are accompanied by a one-loop counterterm diagram, not shown. The curly lines represent gluon 
propagators, and the gluon self-energy bubble shown in diagrams 29-31 includes quark, gluon, and ghost loops. Diagrams 
1-31 contribute to four-quark operator renormalization and are numbered consitently with Refs. [SI IE [|. Diagrams 32-46 are 


considered for the first time here. The 1/e pole structure of each diagram is summarized in Tables III - |VI 
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with equal momenta and gluon exchange between quarks with opposite momenta. With A/ defined by Eq. a 

gluon is exchanged between quarks with equal momenta in 2/5 of the Wick contractions and between quarks with 
opposite momenta in 3/5 of the Wick contractions. Since external quark fields of the same flavor are antisymmetrized 
when constructing A/, each Wick contraction contributes to A/ with equal weight. 

The amplitudes and Aff for diagrams of class d can be evaluated using standard techniques summarized 
in Appendices 0-0 After expressing the resulting spin-color tensors as linear combinations of the basis tensors 
introduced in Appendix]^ it is straightforward to verify that most spin-color tensors contributing to and 
have index exchange symmetries different from the symmetries of the tree-level operator insertion. These contributions 
vanish after making the replacement of Eq. (251 and can be neglected. The remaining contributions to the one-loop 
amplitudes Ai^ for the relevant combination of 2/5 the equal momentum routing amplitude plus 3/5 the opposite 
momentum routing amplitude are given by 




(26a) 




47r \p 


10 


i]r- 


AAS 


(26b) 




1 n 

1^15 60 

+(l0 7/if 0|>7 m) 


+ ^ In 2 + ^ In 2^ ( 7 ^^ 0 |^7 m ® 1 ) 


nASA 


nAAS 


(7m^® 1 ®|>7m) T 


iSAA ^j^AAS 




+ ^ ^ ® ® I'^x^SS (4 0 -k (cr^,. 0 1 0 


(26c) 


1/1 11^ 2 , „ 2k „ 

p2 yi5 60 15 15 




1 [T 


ASA _ '±rpAAS 


) ( ^ rj^AAS 


0 1 0 f' f/x; 

where the color tensors and are defined in Appendix]^ Similarly, the one-loop contributions to Alf with 
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the correct index symmetries are 


/^_^^_4_2^+^ln2+^ln2) [1 0 1 0 


47r \p' 


(27a) 


3 17 2 \ 

-^-^ + gln2j [(ct^^ (g) g i ) t^as ^ ^ g 1 g 


1/1 lie 2 , „ 2 e, „ 




+ (1 0 ® |>7m) ( ) + (7/.^ «) 1 0 |>7m) ( 


+ -T^ 
9 


(27b) 


Mf = 




(-^ - ^ + ^ ln2) [(cr^^ g a^i, g 1)T^'4S ^ ^ + (cr^,. gig (T^i.)T^'^^] 


1 lie 2 , „ 2 e, „ 


in^if ® fi^j. g 1) T 


^AAS 


- -T 


•sss 


+ (1 g g fl^) ( 


(27c) 


Hn^.f ® 10 pi^) ( 


To complete our calculation of ., we need to precisely define operator counterterms that renormalize the vertex 
functions associated with these amplitudes. Subtleties arise at this step. These subtleties and their resolution are 
discussed in the following section. 


B. Evanescent Operators 


In order to precisely define operator counterterms suitable for RI-MOM or MS renormalization, we must address 
the issue that our operator basis is complete in Z? = 4 but incomplete in arbitrary D. This issue also arises for 
four-quark operators in weak matrix element calculations. For four-quark operators it has been consistently resolved 
through the introduction of evanescent operators vanishing in H = 4 |5TJ [51] . Following this approach, in this 
section we precisely define evanescent operator counterterms for the Qi. It would be possible to present complete 
one-loop-matching results without these precise definitions, but the definitions and notation introduced in this section 
will prove essential for calculating the non-trivial evanescent contributions to two-loop running in Sec. jVIj 

Renormalized vertex functions include counterterms that remove the 1/e poles in Eq. (26)-(27). The physical 
operators Qi mix under RG evolution with the operators used to construct these counterterms, so a complete operator 
basis must include all operators used to construct counterterms. Dirac structures involving a ® a appear in the pole 
terms above, but have been eliminated from our complete basis in D = 4 by means of the Fierz transformation 


(sPfpf - AP^P^f) • 


(28) 


This relation follows from completeness of a basis of 16 Dirac matrices in D = 4 and cannot be uniquely continued to 
an analytic function of D. In particular, one could prescribe that in the dimensionally regularized theory 


[CP^a^X^[CP^.a^,V^ ^ ((8 + ai£)P“'P/ - (4 + a2e)PfP^f) . 


(29) 


with ai and 02 arbitrary. The choice oi = 02 = 0 ensures that g is kept equal to its P = 4 Fierz transform. 
Conversely, ® 'Jij.'ju is not kept equal to its P = 4 Fierz transform with ai = 02 = 0. The necessity of breaking 
one or the other Fierz relation follows from the well-known property that contraction of a tensor operator with 
does not commute with renormalization of the dimensionally regularized tensor operator m 
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Working in a dimensionally regularized theory, it is important to distinguish between counterterms in the span of 
the D = A basis operators Qj in D dimensions, and counterterms that are in linearly independent in D dimensions. 
A convenient basis is found by including Qi along with a set of evanescent operators Ej that vanish in £) = 4 but 
are needed as counterterms to renormalize matrix elements of Qi. For example, renormalization of Qi requires a 
counterterm insertion of 


= {'ijjCPRaf,„iT2tp){'ipCPR<Tf,jT2tlj){tpCPRiT2T+'ilj)T^^^ - 12Qi. 


(30) 


D=4 


The evanescent operator Ef vanishes in Z? = 4 by Eq. (281, color-flavor Fierz relations, and the D = A Fierz relation 


Qi = Qi- Because Eq. (28) is broken in dimensional regularization, it is possible that loop-level corrections will 


introduce 0{as) contributions to matrix elements of that do not vanish in H = 4. Explicit calculation demonstrates 
that this possibility is realized. The non-vanishing one-loop contributions are 0(e°) and arise from Xje poles in one- 
loop integrals multiplied by 0(e) suppressed differences in Dirac algebra for the two terms on the RHS of Eq. (301. In 


a naive definition of the MS renormalization scheme, matrix elements of RG evolved physical operators will include 
non-vanishing contributions from renormalized evanescent operators. 

Eollowing Ref. [59], we adopt a definition of the MS renormalization scheme in which the renormalized evanescent 
operators Ej mixing with Qi under RG evolution are defined to include finite 0(e°) counterterms. These counterterms 
are chosen to make loop-level matrix elements of Ej vanish in D = 4 at a particular scale /i. It is proven in Refs. |5TI[5^ 
that this is sufficient to make renormalized matrix elements of generic four-quark evanescent operators vanish at all 
scales. Extension of this proof to six-quark operators is straightforward and discussed in Sec. jVTj The basis used here 
for the Ej needed as one-loop counterterms for Qj is explicitly presented in Appendix Physical observables are 
independent of evanescent basis, but renormalized Wilson coefficients and matrix elements separately are not. It is 
therefore imperative that this same basis is used for loop-level BSM matching calculations. This subtlety is irrelevant 
for tree-level BSM matching calculations. 

Inclusion of MS counterterms in this scheme is equivalent to replacing all terms involving a ® a/e with their D = A 
Fierz transforms plus the evanescent operators of Appendix]^ for instance 

PRa^^iT 2 'il}){'il}CPRa^i,iT 2 i}){'ij^CPRiT 2 Tj^)T^^^ + (MS counterterm) = \l2Qi + z-E'i- (31) 


Inclusion of these counterterms leads to mixing between Qj and Ej and we must enlarge our basis of renormalized 
operators to 


^Qiip)^ 


^ Zii{p) ZjEjip) 



Ei{p) 

= 

Ze,i{p) ZejE,{p) 


A? 

V ^ ) 



■) 

V ^ / 




fQ°\ 

A? 

V ^) 


(32) 


where the ellipses indicate that increasingly many evanescent operators are required to form a complete basis for RG 
evolution at increasingly high loop order. We are specializing to the case of no mixing between the Qj or between the 
Ej, extension to the general case is straightforward. 

The one-loop vertex function Aj^^ can now be expressed in terms of tree-level vertex functions A^^^ built from 

Ej, operator counterterms and vertex functions built from the non-local finite terms in Eq. (26) - (27) 

including and /p?), 




-( 1)4 

^ii 


(1),0 , r(l) 




II 


A 


(0) 


r(l).l 

^lE, 


+ 1 + L\^il'^^i{p, p) + 0{e) 


,(o) 


(1),0. 


(33) 


where all the loop diagram coefficients are pure numbers independent of p, p, and e. and 4>/(p, p) are simply 
obtained by expressing vertex functions constructed from Eqs. (26) - (27) in the form of Eq. (33). The one-loop vertex 
function A^^^ for Ej can similarly be expressed as 


Afi,^ (p, p) 



(34) 


where Fj is a new evanescent operator not included in the Ej that should be included in the ...’s in Eq. (32). MS 
counterterms can be defined as 


_ r(i)4 _ r(i).i _ 

jfMS,(l) _ ^II ;:MS,(1) _ ^lEi ^MS,(1) _ r(l).0 

"// “ g ; °IEi — ^ I ‘^Ejl — ^Ejl ■ 


(35) 
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The bare operator includes UV singularities due to the presence of six bare quark fields as well as the vertex 
function singularities above. In order to remove all UV singularities from (5/(/r), define 




-3 




Itt 


+ 0(Q!s) 


(36) 


where represents two-loop counterterms that will be explicitly constructed in Sec. 

j-MS.(i) / asifJ-) 


VI 


For future use, define 


9.MS 

^lEi 


+O(of) 




MS , 


(a^) = {z, 

(rt=fe" (^)+ c"’ +o(c 


(37) 


/ 7 MS 

^EiEi 


(A^) = {z 


MS 


-3 




.(Ai) 


47r 


+ 


This completes our definition of MS operator renormalization factors in terms of diagrammatic counterterms. 

The RI-MOM operator renormalization condition should also be modified so that RI-MOM renormalized evanescent 
operators have vanishing matrix elements in I? = 4. This is accomplished by adding a supplemental RI-MOM 
renormalization condition 


Ty{ViAej)^0. 


Combining this with the RI-MOM condition Eq. (20) expanded to 0{as) gives 


^RI.(l) ^-^ p = po)] , 


L 


(i).i 

II 


-(i),o 


-(l).Or 


where in analogy to Eq. (36), 


Zf/(p) = (Zf (p)) 


-3 


i+4 'V)(?^)+o(o;) 


(38) 


(39) 


(40) 


The one-loop-matching factor defined in Eq. (231 therefore has the diagrammatic expansion 

^(0) ^ _ ^m(i) _ 3Zfi'(i)(po) + 3ZfS’(i)(po) 

= - L«;°Tr [7^/<i>,(po, A^ = Po)] - 3Zf’(i)(po) + 3Z^’(^)(po). 

Applying the diagrammatic results of the last section gives 

(0)^101 13^ , 8 

31 7^ 88 


1 + 


56$ . 


=-T-¥ + 15'"^ + TV‘''^- 

177 14{_64 

“ 10 + 5 5 

(0) 49 3$ 24 8$ 

5do)__109_13$ 1 8$ 

1 “30 15 + 15 + 3 ’ 

7<3°) = -^-^+^ln2+^ln2. 

3 10 5 5 5 


(41) 

(42a) 

(42b) 

(42c) 

(42d) 

(42e) 

(42f) 

(42g) 


The final one-loop-matching results in Table are obtained after choosing Landau gauge, $ = 0. 
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VI. TWO-LOOP RUNNING 

In order to simultaneously remove large logarithms from perturbative calculations of RI-MOM matching factors and 
BSM Wilson coefficients, RG evolution can be used to relate Wilson coefficients calculated at different renormalization 
scales.^ We perform this RG evolution in the MS scheme for simplicity, and all quantities in this section with suppressed 
renormalization scheme labels are in the MS scheme. The renormalization scale dependence of the Wilson coefficients 
can be determined from the MS anomalous dimension matrix 


luiois) = 


dln/r 




K 


dln/r 


— (0) 
= rjj 


Ltt 


■if) 


dTT 


(43) 


+ 0(ag), 


where the first equality follows from renormalization scale independence of and the second defines the pertur¬ 
bative expansion coefficients and appearing in Eq. (H- The other factors appearing in Eq. ii) are related 
to the QCD /3-function, defined by 


dln/r 


as(/r) = 2j3{as, e)as{y) = (-2e -|- 2/3(as)) Q!s(/r) 


= (_2s-2/3o^-2/3/“^(^) 


(44) 


4 ^ 47 r ' +0(«D)«.(/^)- 

The U-independent piece of the /3-function has a perturbative expansion that is conventionally written as 




4» / 'W 4» ' 

where for QCD with Nf active quark flavors the well-known perturbative coefficients are 


l3o = n--Nf, 


00 

/3i = 102 - -Nf. 


(45) 


(46) 


In MS and other minimal subtraction schemes, fj, dependence of Zjj and y/j only enters through dependence on 
asifJ,). The differential equation in Eq. (43) can be readily solved in a diagonal operator basis where y/j = (J/jy/, 


g/(M2) 

CM 


= exp 


/ 7/(asAi)) — 

/mi ^ 




-7fV(2/3o) 


= exp 


2/3o^ 


y,(a') 

,(mi) 2/3(a() 

Os^ig) - as(Mi) 


2/3o 


47r 


+ C{al) 


(47) 


This equation can be used to RG evolve BSM-scale Wilson coefficients between quark mass thresholds at which the 
number of active flavors Nf decreases. In Sec.|^we introduced as the renormalization scheme matching factor 

relating and Qf^. The renormalization scheme invariance of 'H^Jf allows and to be related using 

rij{p). This allows us to express appearing in Eq. ([^ as 


I CMSr ) CMS( ) 


GrS(^) cfS(/i) 

a.(po) ^ 
as(/i) 




I — r 


An 

( 0 ) ctsjpo) 
An 


Mi 


(0) 


2/3o^ 


2/3o 


o-siPo) - asjp) 
An 


+ 0{al) 


( 48 ) 


^ See Ref. | 74| for a nice review of RG evolution for weak matrix elements including discussion of evanescent operators. 
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The remainder of this section discusses the diagrammatic evaluation of 7 /^^^ ■ 

In Sec. |VB| we discussed the need to remove dimensional regularization artifacts by adding finite counterterms 
proportional to evanescent operators Ej to physical operators Qj and vice versa. Without these counterterms the 
renormalized Ej would contribute to physical observables and therefore to BSM matching calculations of Wilson 
coefficients. With these counterterms, Qi mixes under renormalization with Ej and the assumption of a diago¬ 
nal anomalous dimension matrix taken above is invalidated. We must instead consider the renormalization scale 
dependence of the infinite dimensional matrix of Eq. (32) and define 

^ 7 /J lIEi ■ ■ ^ 

7 = lEiI lEiEi 

V ^ ■7 


(49) 


Eq. (471 is preserved if and only if ^e,i = 0 to two-loop order. A proof that 7 ^ 7 / vanishes to all-orders for generic 

and applies to our six-quark operators as well. This is discussed in 


four-quark operators is given in Refs, 
detail at the end of this section. 

Since /i dependence of Z only enters through explicit dependence on as(/i), the anomalous dimension matrix 7 is 
given by 


dn 


7 = - ( ] ■ Z = - ( 2l3{as,e)as{fJ,)-^Z ) • Z 


d 


7-1 


(50) 


where • denotes matrix multiplication. Perturbative coefficients of Z defined in analogy to Eq. (43) are given by 

7 ^°) = 2eZ^^\ 

7 ( 1 ) = 4e7(2) - -h 2l3oZ^^l 

The anomalous dimensions of the physical operators Qi are therefore 

7^^ =2 ezW, 


(51) 


7i^^ = 4£Zf; - 2s (Zjf + Z^^Z^j^l) + 2p,zf}. 


(52) 


The non-trivial effect of evanescent counterterm subtraction is the appearance of in 77^- 

The factors above are simply related to diagrammatic counterterms. The one-loop anomalous dimension is deter¬ 
mined by the counterterms of Sec. lYl as 


( 0 ) 

7 ) 


= 2s {dP - 3^7)) 


(53) 


which is finite at I? = 4. Calculation of 7 ^^ requires the 1/s pole contributions to the two-loop Qi vertex functions 




72).2 

^II 


d 2 ).l 

^II 


Af 


.( 2 ), 2 ' 

^lE, 

e 2 


A' 


(0) 

E, 


.( 2 ). 2 ' 

^IFi 

e 2 


A' 


(0) 


0 (e° 


(54) 


Including one-loop counterterm diagrams with insertions of 6^ as well as quark self-energy, gluon self-energy, and 

( 2 ) 1 

quark-gluon-vertex counterterms ensures that non-local divergences are cancelled and A//’ is a pure number. The 
two-loop MS counterterm is then defined as 


7(2).2 

d 2 ) _ _7li_ 

“ 7=2 


72).1 
^II 


(55) 


We can then use Eqs. (36) - (37) to express the Z factors appearing in Eq. (52) in terms of these counterterms. 


71) = 4e (47 - 3^7)) - 2s (d 


( 1)2 

II 


4444 - 3^4^') + 2/3o (47 - 37'0 ■ 


(56) 


Two-loop counterterms include l/s^ contributions, so the various terms in Eq. (56) are divergent in A) = 4. Renormal- 
izability of composite operators in the MS scheme guarantees that matrix elements of Qi{fJ-) are free of UV divergences 





















18 


at all renormalization scales and therefore that 7 / is finite order-by-order 177]. This means that divergences must 
cancel between the terms of Eq. (56).® After this cancellation, the anomalous dimension is given by 


qf) = - 12Zf ) 


(57) 


where is the 1/e piece of Eq. (16). 

It was noticed in Ref. |59j that the finite contributions to 7 }^^ from mixing with evanescent operators contribute ex¬ 
actly like an additional counterterm diagrams apart from the relative factor of (— 1 / 2 ) between and in 

Eq. (57). As discussed after Eq. (54), includes contributions with one-loop counterterm diagrams containing inser¬ 
tions of S^j^j . Suppose for each of these one-loop counterterm diagrams we include an additional counterterm diagram 
with an insertion of (l/2)i5^^^j = (—1/2)L^^^}°. These diagrams make a 1/e pole contribution of (—l/2)Lj^(^L^^^}°/e. 
Including these additional one-loop counterterm diagrams with insertion of (1/2)^^^^^ therefore shifts the 1/e single¬ 
pole part of to 


rr 1 ( 2).1 


= L 


( 2),1 

II 


_ 1 r(l)4;-(l).0 
2 lEi ^Ejl ’ 


(58) 


the factor appearing directly in Eq. (57). 


To ensure a proper treatment of two-loop subdivergences and verify cancellation of non-local divergences diagram- 
by-diagram, each two-loop diagram should be combined with a one-loop counterterm diagram in which any divergent 
one-loop subdiagram present is replaced by a one-loop counterterm that cancels the subdivergence. To provide 
this cancellation diagram-by-diagram, the one-loop counterterm must have the same color structure as the one-loop 
subdiagram. For subdiagrams with the topology of a one-loop self-energy or vertex correction, the color structure of the 
subdiagram is a simple multiple of the corresponding tree-level color structure, and this procedure is straightforward. 
For subdiagrams with the topology of a one-loop operator correction d = 1-3, the color structure of the subdiagram 
differs from the tree-level operator color structure, and care must be taken. In particular, the physical operator 
counterterm associated with each d = 1-3 subdiagram must precisely reproduce that subdiagram’s contribution to 
in D dimensions. This means that the operator used for each physical operator counterterm must be proportional 
to Qi in D dimensions. 

To produce physical operator counterterms proportional to Qi in D dimensions, it is necessary but not sufficient 
that all appearances of cr 0 cr in d = 1-3 subdiagrams are replaced in one-loop counterterm diagrams by trivial spin 
structures containing only the identity matrix. A convenient prescription that meets this necessary criterion is to use 
Eq. (29) with ai = 02 = 0. This prescription amounts to replacing all appearances of cr 0 cr in d = 1-3 subdiagrams 


with their D = A Fierz transforms and clearly produces one-loop counterterm diagrams with the same color structures 
as each two-loop diagram. However, one-loop counterterms defined by this prescription are not proportional to Qi 
in D dimensions and instead differ by terms proportional to (Qi — Qi) and (Qa — Qs)- Applying the prescription of 
Eq. (29) with ai = 02 = 0 to the one-loop amplitude defines a basis E/ different from the one-loop counterterm basis 
E[. For example, the operator A/' is given by 


{'ii)CPRa^^iT2'i/){i/CPRa^uiT2'/i){'/iCPRiT2T+ti))T^^^ 

= i{^/^[CPRrhT2ilp{r[CPRV^ir2p){.r[CPRplT2T+^^)T^^^ 
= 12Qi — 6 (Qi — Qi) + E^', 




(59) 


where the last equality follows from color-flavor algebra. Operators E/', E^', E[, and Eg differ from their unprimed 
counterparts by factors of (Qi —Qi) and {Q3 — Q3) as seen for A/' by comparing Eq. (59) and Eq. (30). The remaining 
operators needed to define the E/ basis and it relation to the Ej basis are explicitly presented in Appendix O 

In the E'j basis, the total physical plus evanescent one-loop operator counterterm associated with two-loop diagrams 
d = 4-6, 16-24, and 32-45 containing d = 1-3 subdiagrams is simply given by minus the pole part of the subdiagram 


A potential point of confusion: if one naively takes Eq. {5TJ with Z replaced by Zg ^ as a formula for the two-loop quark field anomalous 
dimension and inserts Eq. ( |16| |, 1/e^ divergences do not cancel. The subtlety is that Zq depends on the gauge parameter which in 
turn depends on the renormalization scale Ezj. When Eq. l |50[ > is modified to include this additional source of renormalization scale 
dependence, the resulting quark field anomalous dimension is indeed finite. 
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with appearances of tr ® tr replaced by 




I a + F[a 0 cr]) 


(60) 


where F[a ® a] is given by the RHS of Eq. (29) with oi = 02 = 0. The total 1/e single-pole contribution given by a 

(2') 1 1 

diagrammatic two-loop calculation using this prescription we denote and is related to [Ltot]ii’ by 

1 


rr 1 ( 2 ),1) 


= Wtot] 


(2),1 

II 


-r(l) (r(l).O .(1).0 

2^IEi \^E,I 


(61) 


The change of evanescent basis factors appearing in Eq. ( |61[ ) can be immediately obtained from one-loop results for 
~^*^^) E^nd and are given in Appendix]^ We have explicitly verified that after including these change 

of evanescent basis factors the total contribution of one-loop physical operator counterterm diagrams is equal to the 
sum of one-loop counterterm diagrams with insertions of Sj^j . 

With this evanescent-counterterm diagram prescriptions in hand, diagrammatic calculation of the E/ basis contri¬ 
butions to 7 ^^^ proceeds as described in Sec. V A and Appendix | a]-|d| There are 320 contributing two-loop diagrams 
organized into independent classes d = 4-46 in Fig. Diagrams 4-31 contribute to NNLO renormalization of scalar 
four-quark operators, and many also contribute to renormalization of three-quark operators. We have adopted the 
same numbering scheme for these diagrams used in Ref. [Ml inn, and have verified that our results for these diagrams 
agree with the scalar four-quark operator results of Ref. m after changing to the appropriate evanescent operator 
basis (the “Greek projection” basis, see Appendix]^. Diagrams 32-46 are new. For future calculations, it is interesting 
to note that there are no additional two-loop diagram classes appearing for operators with more than six quarks. In 
principle the two-loop anomalous dimension of any AB = N operator composed of a product of scalar diquarks could 
be computed from the results of Table m combinatorics, and group theory. 


After including all appropriate one-loop counterterms, including the E/ evanescent counterterms defined by Eq. (601, 


the 1 /e pole part of each diagram can be decomposed into a color factor for each diagram times a linear combination 
of Dirac structures that is identical up to diquark permutations for all diagrams in the class. The pole parts of these 
Dirac structures are shown in Table |III1 and the corresponding color factor for a representative diagram is shown 


in Table IV The 1/e pole coefficients of the combined spin-color tensors with non-vanishing contributions to 
operators are shown in Table The corresponding pole coefficients contributing to T'®'®'® operators are shown in 


Table VI After re-introducing quark fields and flavor tensors for a given operator by Eq. (25), the resulting pole 


structures for each operator are related by Eq. (Cl) to a multiple of the original operator plus irrelevant contributions 

When the dust settles, these diagrammatic contributions sum to L^^j’^/e"^ + We have explicitly verified 

that the 1/e^ contributions to cancel with the other divergent terms in Eq. ([5^. This provides a highly non-trivial 


check on the calculation. The physical anomalous dimensions are then given by Eq. (57) as 


( 1 ) 335 34V/ 

“ 3 9 ’ 

(62a) 

( 1 ) _ 91 26V/ 

3 9 ’ 

(62b) 


(62c) 

,<■>- 229 

(62d) 

7 W ^ 238 - 14V/, 

(62e) 

^( 1 ) 797 1I8V/ 

“ 3 9 ’ 

(62f) 

7C’-218D8«/^ 

(62g) 


In the MS scheme used in this work, Fierz-conjugate operators Qi , Qi and Q 3 , Q 3 are equal in D = 4 but do not have 
identical two-loop anomalous dimensions. Conversely, to be regularization independent the RI-MOM scheme cannot 
assign different anomalous dimensions to operators identical in D = 4.® Two-loop RI-MOM anomalous dimensions 


^ We thank Sergey Syritsyn for bringing this point to our attention. 
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are gauge dependent and depend on the external state appearing in the RI-MOM renormalization condition, but 
g] 2 ould be independent of evanescent basis and equal for Fierz-conjugate operators. The two-loop RI-MOM 
anomalous dimension 


(l).RI (1),MS , oo (0) 

7 / =7/ +2/3or}^ (bdj 

can be shown to be independent of the renormalization scheme and evanescent basis used to define counterterms |741 
[HZ]. This provides a pair of consistency conditions 


RI,(l) ~RI,(1) RI,(1) ~RI.(1) 

7i =7i , 73 = 73 > 


(64) 


that can be readily verified to hold for the results of Eq. (42) and Eq. (62) in a general gauge. Eq. (64) provides 
a useful check on our calculation, and particularly on evanescent contributions to that cannot spoil the 1 /e^ 
cancellation consistency check. In particular, we observe that violations of Eq. (641 and 1 /e^ poles cancel independently 
from the following sets of diagrams: {4-6, 16-24, 32-45} containing d = 1-3 subdiagrams; {7-9} containing crossed 
gluon lines; {28,46} containing three-gluon vertices and no divergent subdiagrams; and {10-15, 25-27, 29-31} combined 
with the one-loop contribution 2 fiQr^p and one- and two-loop wavefunction renormalization. 

It remains to verify that ^EiI = 0. At one-loop order 




1e,i 


Ejl- 


Since = 


is a finite counterterm, trivially vanishes at D = 4. At two-loop order 

IeI = 4^^ - 2^ + S'eU/eI) + 2/3o4V 


(65) 


( 66 ) 


Vanishing of Je'Ii trivial. Refs. [8TJ|82] prove that for generic four-quark operators a strictly stronger statement 

is in fact true: the analog of 7 is upper-triangular to all orders. The argument of Herrlich and Nierste in Ref. [Hlj is 
quite general and only relies on cancellation of non-local divergences and 0 {e) suppression of diagrams with evanescent 
operator insertions. We briefly review their argument to demonstrate that it applies to our six-quark operators without 
modification. 

The evanescent operators Ej vanish in I? = 4 by spin-color-flavor Fierz transformations. The one-loop counterterm 
therefore includes 0(e)-suppressed tensor algebra and can only include contributions non-vanishing in = 4 from 
terms that have received a 1/e enhancement from loop integrals. is therefore O(s^), and one-loop counterterm 
diagrams with insertions of only make 1 /e pole contributions to JejI from terms that recieve additional 1 /e 
integral enhancements. This implies that the only 1/e pole contributions to 7 ^^^^^ from and S^eIei^e^iI single 

poles and arise from integral contributions with 1 /e^ double-pole enhancements and 0 (e) tensor-algebra suppression. 
The sum of these one-loop counterterm diagram contributions is explicitly given by /p'^Y(S^ eIjS^'^ + ^e^jEi^e^iI ~ 
^E^ji^o/s)- Without the 0(e)-suppressed tensor algebra, this expression would include non-local divergences arising 


from 1/e^ factors multiplied by /p^Y. Cancellation of non-local divergences is independent of the tensor structure 
of operator insertions, and these would-be non-local divergences must cancel with 1 /e^ two-loop integrals multiplied by 
{p?/p^Y^ and the same 0(e)-suppressed tensor algebra factors. This gives ~^^E^jE,^E^ji~^Ehl^o/£), 

which when inserted in Eq. ( 66 ) gives = 0. The remaining inductive step needed to prove that 7 is upper- 
triangular to all orders does not rely on a particular evanescent operator definition jSTl [82] and applies here as well. 


VII. PHENOMENOLOGICAL APPLICATIONS: AN ILLUSTRATIVE EXAMPLE 

The phenomenological consequences of neutron-antineutron operator renormalization are encoded in the effective 
Hamiltonian R/ff of Eq. (U and the operator renormalization factors 7 |°\ 7/^\ collected in Table I These 

operator renormalization factors govern the relations between matrix elements of Qj with different renormmization 
scheme and scale choices. Non-perturbative lattice QCD determinations of the renormalized QCD matrix elements^° 


10 


Lattice matching scales of po — 2 GeV are typically large enough for matching to be perturbative but small enough that unphysical UV 
cuttoff effects are minimal {Aqcd < Po < where a is the lattice spacing). 
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( 


Qi 


n) can be combined with these operator renormalization factors to determine QCD matrix elements 
at high scales /i, where BSM physics is usually assumed to be perturbative. Once high-scale QCD 
matrix elements have been calculated with fully quantified uncertainties, perturbative BSM matching for a particular 
BSM theory of interest can be used to predict nn transition matrix elements 


^_^Po) 


- = Sm = (n| y-eff , 

'^nri 


(67) 


in terms of basic BSM parameters. Experimental constrains on the neutron-antineutron vacuum transition proba¬ 
bility Pnn{t) = sin^(|(5m|t) then unambiguously constrain the parameter space of BSM theories predicting neutron- 
antineutron transitions. 

The extraction of phenomenological predictions from Eq. ([^ is best explained via an example of one specific BSM 
model. Several broad classes of simplified models with classical baryon number violation but no proton decay were 
recently discussed by Arnold, Fornal, and Wise |15j. The BSM field content of these models consists of a pair of colored 
scalar fields that carry non-integer baryon or lepton number. For illustrative purposes, we will use the simplified model 
discussed most heavily (Model 1 in Ref. [5S]), which we’ll call the AFWl model, to perform our calculation. Identical 
steps can be used to make predictions for other, more complicated BSM theories once a six-quark effective Hamiltonian 
■Hgjy has been determined for those theories. 

The AFWl model adds two new scalars to the standard model, Xi and X 2 , which transform as Xi G (6,1, —1/3) 
and X 2 € (6,1,2/3) under SU{3)c x SU{2)l x U{1)y- The Xi and X 2 couplings to the SM right-handed fermions 
are given by g[ and g 2 , respectively, and an additional three-scalar coupling between two-Xi and one-Ai 2 is given 
by A. This model allows neutron-antineutron transitions at tree-level. The Hamiltonian operator is found by 

evaluating a tree-level Feynman diagram connecting six external quarks all carrying zero momentum. The resulting 
P-eff for the AFWl model is presented in terms of fixed-flavor quark fields uf, df in Eq. (12) of Ref. [25], neglecting 
scalar couplings to left-handed quarks for simplicity. In the fixed-flavor and chiral operator bases, this Hamiltonian 
is given at tree-level by 


4M/M| 


Qi + 

5 


( 68 ) 


where g'l^ = g\^ are dimensionless couplings assumed to be 0(1) at a high scale M and A, Mi, M 2 massive couplings 
assumed to be 0{M). Perturbative corrections to this expression include h\{g?/M"^) factors, and so to allow the 
validity of tree-level BSM matching just described fj, = M is chosen. RG evolution is simplest in minimal subtraction 
schemes such as MS and, as a result, we formally prescribe that these corrections should be calculated in the MS 
scheme. With these renormalization choices and BSM naturalness assumptions, the AFWl Hamiltonian can be 
expressed as 


njnn _ ^ 

~ 16M5 


Qf®(M) + 

5 


(69) 


The operator renormalization results of this work then allow the AFWl Hamiltonian to be expressed as 


njnn _ 

“ 16M5 


C/4(M,po)Qf (Po) + ^UiiM,po)Qf\po) 


(70) 


where 17/(^,po) is the RG evolution and renormalization-scheme-matching factor appearing in Eq. ([^. Explicit 
evaluation of Ui{p,,po) for arbitrary p, po requires an accurate parameterization of as(p). For this we take the four- 
loop parametrization of as{p) in terms of A)lJ® and known /3-function coefficients presented in Ref. |75|. The full RG 
evolution between p and po is included through a product of factors C/^-^ {pi, P 2 ) where Nf is varied across each quark 
mass threshold. Implicit Nf dependence in the parametrization of as{p) in terms of the fit parameters must be 

included along with explicit Nf dependence in /3o, /3i, and 

Preliminary lattice QGD neutron-antineutron simulations have been performed |54] . Updated (but still preliminary) 
results on the two operators of interest are given by |5S| 


(n| Qf\po) \n) = (0.00 ± 2.06) x 10"^ GeV®, 


(n 


Q?\po) 


= (-56.13 ±2.42) X 10"® GeV®, 


(71) 


where the errors shown are purely statistical and fitting errors. It is important to note that these preliminary results 
include several significant sources of systematic uncertainty that have not been quantified at this time. These sources 
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of systematics are the absence of RI-MOM non-perturbative renormalization,^^ unphysically large quark masses (pion 
masses of roughly 400 MeV), lattice spacing artifacts (which unphysically break chiral symmetry), and finite spatial 
extent artifacts. All of these systematic uncertainties can be quantified and reduced with increased computing.For 
this reason, these results should only be viewed as an illustrative example of how to combine the perturbative QCD 
operator renormalization results of this paper, BSM calculations of Wilson coefficients, and non-perturbatively renor¬ 
malized lattice QCD matrix elements to arrive at a physical quantity that can be measured/bounded by experiment. 

The experimental limit on 5m determined from Super K measurements of tqw is Ems] 

\dm\ < 2 X 10"^^ GeV. (72) 

Ref. |25] uses this limit and an estimate for the QCD matrix elements to relate this to a limit on the AFWl model 
scale, 


M > 500 TeV. 


(73) 


Using the operator renormalization factors of Table |T] and Eq. ( [tT] ), it is possible to express matrix elements of Eq. (70) 
in terms of M and known parameters. Substituting M = 500 TeV into Eq. (70), 


\5m\ = (^ -2.44 -tO.33 -O.M 

LO NLO NNLO matching NNLO running 

= (3.68 ±0.16) X lO-^-^ GeV, 


±1^ ) X 10 

Lattice Statistical 


-34 


GeV 


(74) 


gives rise to a nn vacuum transition time for the AEWl model more than five times longer than the Tnn predicted 
by the QCD matrix element estimate used in Ref. [35] . Note that NLO one-loop running provides a multiplicative 
correction to the LO matrix element and is shown as an additive correction above only to illustrate the size of the 
correction for this example. Corrections beyond NLO can be organized a perturbative power series in as{po). This 
perturbative series for \5m\ appears to be converging nicely, with O{as{po)) NNLO corrections changing the NLO 
result by 7%. Assuming further corrections have the same rate of convergence, we expect unknown O{as{po)'^) N^LO 
corrections to modify the NNLO result by ~ 2%. 


A constraint on M can be derived by inverting the experimental bound Eq. (72) 


M>(^ -34 ±5 -16 

LO NLO NNLO matching NNLO running 

M > 357 zb 3 TeV. 


^ ) TeV, 

Lattice Statistical 


(75) 


This constraint is nearly one-third weaker than the constraint estimated in Ref. |25j . 


VIII. CONCLUSION 

To determine which BSM theories are able to produce the observed baryon asymmetry of our universe, it is essential 
that each make reliable predictions for CP violating and baryon-number violating processes that can be probed 
experimentally. Theories with AR = 2 interactions can provide viable baryogenesis mechanisms while avoiding 
stringent experimental bounds on AR = 1 proton decay rates. Some of these theories predict new physics in the 
100-1000 TeV range that induce nn transitions that are just outside the reach of current experimental bounds. If 
these theories can be reliably constrained by experimental measurements of the nn vacuum transition time t„h, next- 
generation nn experiments can search for new physics appearing at scales comparable to, or higher than, scales probed 
in next-generation collider experiments. 

Reliable predictions for 5m = \lTnn-, the parameter governing the neutron-antineutron vacuum transition probabil¬ 
ity Pnn = sin^(|(5m|t), can be made by perturbatively matching BSM theories to an effective field theory containing 
Standard Model operators. The six-quark operator matrix elements contributing to 5m can be calculated with lattice 
QCD at computationally accessible lattice matching scales of 2 GeV, and large logarithmic strong interaction cor¬ 
rections can be included using perturbative operator renormalization. At NLO, operator renormalization introduces 


Renormalization approximated by tapdole improved tree-level renormalization 18911901 . 

The required computational resources are expected to be similar to those used for lattice calculations of Bk with physical pions and 
chiral fermion discretization m- 
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known multiplicative corrections to six-quark operator matrix elements m- Further operator renormalization cor¬ 
rections are organized as a perturbative series in which the largest contributions arise from NNLO two-loop-running 
and one-loop-matching effects. These effects are calculated for the first time here and summarized in Table |Tj In addi¬ 
tion, operator projectors needed for non-perturbative renormalization of nfi operators and the chiral transformation 
proper ties of nfi operators are presented. 

Sec. VII discusses the calculation of 5m in a simplified model from Ref. [25] in order to illustrate how perturbative 
operator renormalization results are combined with lattice QCD six-quark operator matrix elements and experimental 
bounds on 5m to constrain the scale of new baryon-number violating physics. For fixed BSM parameters, the nfi 
vacuum transition time calculated with the perturbative operator renormalization results of this work and preliminary 
lattice QCD results is found to be more than a factor of five longer than was previously estimated. Several features 
of this simple example have generic implications for more complicated BSM theories and deserve explicit mention: 


• is described by a linear combination of multiple chiral basis operators that make very different contributions 
to nfi matrix elements. Color and flavor structure matters when calculating the nfi vacuum transition time 
predicted by a particular model. 

• (54 has a large positive anomalous dimension and its contributions to 5m are suppressed compared to other 
operators. 

• Q 2 is the only chiral basis operator with a negative anomalous dimension and its contributions are enhanced 
compared to those of other operators. It does not contribute in the simplified model considered. 

• Qi arises in real models. Qi and Qi should be treated on equal footing, and (5 i^^(m) 7 ^ must be 

remembered during BSM matching calculations. Identical considerations apply to Q 3 and Q 3 . 

• Qsj Qsj s-nd Q 7 violate electroweak gauge invariance. They do not appear in the simplified model considered. 

• Assuming can be expressed as linear combinations of Qi,..., Qy, Qi, Q 3 and their parity conjugates 

without the use of (spin) Fierz relations, evanescent operators do not need to be explicitly included in tree-level 
BSM matching calculations. 

• At a fixed BSM scale of 500 TeV, NNLO effects correct the NLO-I-LO 5m prediction by 14%, within the generic 
range of < 26%. Assuming that further perturbative corrections have the same rate of convergence, unknown 
N^LO corrections are estimated to change the NNLO-I-NLO-I-LO 5m prediction by 2%, within the generic 
estimate of < 7%. 


• Perturbative corrections to BSM scale constraints are smaller than corrections to 5m. In the simplified model 
considered, NNLO effects change the NLO-I-LO BSM scale constraint by 3%. 


Without knowledge of NNLO two-loop-running and one-loop-matching factors, perturbative operator renormaliza¬ 
tion effects contribute large unquantified systematic uncertainties to BSM nfi vacuum transition time predictions. 
Including NNLO effects and estimating the size of unknown N^LO corrections turns perturbative operator renormal¬ 
ization into a few-percent-level uncertainty. This places perturbative QCD corrections to Tnn firmly under control. 
Electromagnetic one-loop-running corrections have also been calculated in Ref. m. though a complete electroweak 
one-loop-running calculation has not be performed. 

A complete lattice QCD determination of the nfi matrix elements with controlled systematic uncertainties is nec¬ 
essary to remove the largest remaining Standard Model uncertainties present in BSM predictions of Tnn [MIIHH]- In 
particular, RI-MOM scheme non-perturbative renormalization factors should be calculated, continuum and infinite 
volume extrapolations should be performed, and nfi matrix element calculations should be repeated with physical or 
near-physical pion masses. All of these systematic uncertainties can be removed using existing lattice QCD technology 
and computational resources available in the near future. 
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Appendix A: Tensor Algebra 


This appendix presents Fierz-type relations useful for resolving complicated spin, color, and flavor tensors in a 
desired tensor basis. All relations are derived using a well-known tensor reduction strategy: write the tensor under 
consideration as a linear combination of chosen basis tensors .. .with unknown coefficients ci, C 2 , ..., e.g. 

= c^Bf + C2Bf. (Al) 


It is often useful to chose basis tensors with definite index exchange symmetries. Contracting both sides of this 
equation with each basis tensor gives a system of equations 


/ / BfBf BfBf\ (Ci\ 

^BfBf BfBfj yczj ’ 


(A2) 


that can be readily solved for ci, C 2 . 


1. Color Algebra 


There are five independent su(3)c' tensors that can combine six quarks into color singlet operators, 


Tf 


Tr" 


^ [ij] [kl]{mn} ’ 

These hve basis tensors are constructed from 

-^sss 


T, 


ij]{kl}[mn] 


T> 


{ij}[kl] [mn] 


rpA.. 

^ [ij][kl][mn]- 


'^{ijBkl}{‘mn} — ^ikm^jln “t“ ^ jkTn^iln “t“ Silm^jkn “t“ ^ikn^jlmt 


T, 


AAS 

ij\ [kl] {mn} 
nAAA 


— ^ijm^kln H” ^ijn^klmi 


(A3) 


(A4) 


'^[ij][kl][mn] ~ ^ijm^kln SijnSklm, 


where Sijk is the completely antisymmetric Levi-Civita tensor for su(3). Each tensor is symmetrized { } or anti¬ 
symmetrized [ ] in the three index pairs shown. This corresponds to combining the six 3 quarks as products of 
symmetrized 6 or antisymmetrized 3 diquarks. These tensors also obey the diquark exchange symmetries 


fj^SSS 


= Ti 


<sss 


{ij}{kl}{mn} {kl}{ij}{mn} 


= T\ 


•sss 




T, 


AAS 

ij] [kl]{mn} 


= T, 


AAS 


[kl] [ij]{mn} 


rpA 


ij] [kl] [mn] 


= -T, 


[kl][ij][mn] 


= -T, 


AAA 

ij] [mn] [kl]' 


The remaining basis tensors defined by 


rpASA _ rjiAAS 

[ij]{kl'\[mn] [ij][mn]{kl'\ "> 


rjiSAA _ rpAAS 

^ {ij}lkl][mn] ~ ^ [mn][kl]{ij}' 


(A5) 


(A6) 


When evaluating Feynman diagrams, one encounters contractions of the color tensors and present in 

Qi with the su(3) generators The resulting color tensors can always be expressed in terms of index permutations 
of the original color tensors. For most diagrams this is accomplished through the textbook identity 


* * j j 


= X ( Si'jSj'i — 




(A7) 


where we assume the normalization Tr(f'^t^) = . Certain classes of diagrams involving three-gluon interactions 

require the additional identity 


cABC-lA ,b ,c 
^i'i^j'j ^k’k 




\^k’j ^i'j^j'k^k'i) 


(AS) 


which we have derived by performing the tensor reduction of Eq. (A2) for a basis of Kronecker-delta products. The 


color structure produced by any diagram can therefore be determined from the relations of the generators above and 
color Fierz identities relating index-permuted tensors to the hve basis tensors. 
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The symmetrized color tensor obeys the Fierz identity 


rpSSS _ _^rpSSS _ ^rpAAS 

^ {kj}{il}{mn} 2 ^ 2 ^ [b1 [fc*l{"i"-} ' 


(A9) 


The corresponding relations for interchange of any other index pair follow from the symmetries above. The mixed 
symmetry color tensor obeys the Fierz identities 

rjiAAS _ ^rpSSS I ^rpAAS 

rpAAS _ 1 rj^AAS 1 rpASA I rpAAA 


All other index exchange relations follow by symmetry. The antisymmetrized color tensor obeys the Fierz identity 


rpAAA _ ^rpASA ^ rpSAA 

-‘-[kj][ir\[mn] ~ 2 ^ 2 ^ ’ 


(AlO) 


with all other relations again following by symmetry. 

Color factors produced by any diagram can be expressed in this basis through repeated application of these Fierz 
identities or alternatively by direct tensor reduction of color factors in the forms of Table |IV[ We hnd the second 
approach more convenient at the two-loop level because it can be readily automated in a computer algebra program 
such as Mathematica. 


2. Dirac Algebra 


QCD loop diagrams introduce additional factors of • • • into each diquark. In a theory with massless quarks 
perturbative corrections will not modify the chirality of each diquark. We can therefore express any Dirac structure 
produced by a loop diagram as CP^J^i ® CP^^T 2 ® CP^^T^, where the chirality labels are identical to those of the 
tree-level operator in question. In D = 4, a suitable basis of chirality preserving ri®r 2 ®r 3 independent of quark 
momenta is given by 


I (g) 1 (g) I, ® (g 1, I (g (g cr^jy, cr^jy (g 1 (g *cr^jy (g cr^p (g cr,yp, (AH) 

where = §[7p,7;y] and 1 represents the 4x4 identity matrix. An additional independent structure (g 

Cpp ® CTjyT is not produced in the diagrams considered here. When discussing these basis tensors we will often omit 
the Lorentz indices and write shorthand expressions like cr (g cr (g 1 and a ® a ® a. 

These basis structures provide a convenient orthogonal basis for tensor decompositions of two-loop Dirac structures. 
Operators built from these basis structures are not explicitly included in our physical operator basis. Using spin Fierz 
relations, each basis tensor can be related to a combination of index permutations of 1 (g 1 g I and therefore to the 
physical basis structures. Different techniques are required to find basis decompositions for Dirac structures produced 
in loop diagrams that are valid in D dimensions. Useful discussions on the Z?-dimensional Dirac algebra needed for 
three- and four-quark operator renormalization can be found in Refs. [671 [721 [741193], and in particular many results 
and techniques for D-dimensional tensor reduction can be found in Refs. [HIIIM]. 

As discussed at length in Sec. |VB[ spin Fierz relations are broken in dimensional regularization since the 16 Dirac 
matrices {1, 75 , 7 p, CTpiy} (with ^ < v) are not a complete basis for Dirac matrices in general D. Spin Fierz relations 
should instead be considered prescriptions for defining evanescent operators built from the difference between the LHS 
and RHS of these identities, see Refs. [51171]. One has the freedom to add 0(e) terms in defining this prescription, 
for example 

[CP^apX^[CP^,apA'^^ ^ ((8 + aie)[OP^]“'[OP^]^^ - (4 + a 2 e)[OP^]“^[OP^]'^') . (AI2) 


The 0(e°) coefficients can be calculated by performing a tensor reduction in P = 4. Our basis of evanescent operators 
is explicitly defined in Appendix [^ When applying prescriptions such as Eq. (AI2) to define evanescent operators, 
finite matching factors and Wilson coefficients will depend on the chosen 01 , 02 . Basis dependence cancels so that 
physical quantities such as P"// independent of oi, 02 . Relations between renormalized matrix elements calculated 
with different one-loop evanescent bases follow from general considerations of renormalization scheme dependence, 
see Refs. [60l|82]. Alternative schemes for defining evanescent operators can be found in Refs. [67l|72l|8l]. 

Additional Dirac structures appear in two-loop diagrams that are independent in general D. These must be treated 
with analogous spin Fierz evanescent operator prescriptions, such as 


CPx ^ CCT^i/<7p-Y y ((48 -k bie)CP^ ® CP^ -k (8 -k h2e)CP^ ®CPx ^py ) ; 


(AI3) 
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where bi and 62 are arbitrary parameters used to specify a basis for two-loop evanescent counterterms. Freedom to 
specify 61,62 and other two-loop spin Fierz prescriptions suggests there is an additional ambiguity in 7 ^^^ besides 
the choice of 01,02 that determines the one-loop evanescent counterterms. This suggestion is false.Since one- 
loop-matching factors are independent of the 6 ’s, there is no way for 7 ^^) to depend on the 6 ’s while keeping 

independent of this arbitrary basis choice. Independence of 7 }^^ on the 6 ’s is proven for four-quark operators in 
Ref. [12] • The proof only relies on cancellation of non-local divergences and the factor of 1/2 multiplying evanescent 
counterterm diagrams in Eq. (601 and applies to our six-quark operators without modification. We have explicitly 


verified that cancellation of 6 ’s dependence occurs diagram-by-diagram between two-loop diagrams and one-loop 
evanescent counterterm diagrams in our calculation. 

In addition to Eq. (A13), two-loop diagram evaluation requires the D = A spin Fierz identities 


CPx. <ypT ^ py (4CPx «) CP^ ^ py ' 

-AiCP^a^y I 


.CP 




.CP^ 

^ n n ®CP^ 

^yp) 5 


(A14) 


where 


^X — '^XlX2'^X2X3 


(A15) 


vanishes unless all three diquarks have identical chirality. Relating the above Dirac basis tensors to permutations of 
1 (g) 1 (g) 1 also requires 


^[CP^,a,.]“^[CP,,a^,]'^^[CP, 3 a.p]''^ (l 6 [CP^]“«[CP^]^^[CPJ''^ 

- 8[CP^]“'[CPj^^[CP^]’'f - 8[CP^]“''[CP^]^«[CP^]’'' 
- 8 [CP^]“«[CP^]'^^[CP^]^'' + 8 [CP^]“^[CP^]^'[CPJ’'«) 

Other useful formulae are derived by combining Eq. (A12) - ( A16| ). A particularly useful identity is 

- al3-y6r}(^ 




1 . 1^1 

-l(g)cr(g)cr-|--cr(g)l(g)cr — -crOcrOcr 


Fierz relations involving ^ are also useful when computing evanescent counterterm diagrams 

^XX'® {P'Pxfltj) ~ ^XX'i^Px) ® i^Px) + Px'^ ® (C'PxO’^ii/), 
^{CP^a^^^xfj (g) {CP^^a^^f^x) °=\25^^>{CP^) O (CP^O - (CP^a^^) O {CP^^a^y), 
^{CP^a^y"fxf) O {CPy^'fyxC7f_,y) ^=h2Sx^>{CP^) O (CP^O -I- i{CP^a^y) g) {CP^'a^y). 
Additional Fierz relations are useful for diagram classes 34-45, in particular 


=4l 


1 , 


Pf[Px'i,tV'[Pxfi,T^ =\[Pxi,tr^[Px'ih^^'P^^ + -APx^^^uipfr^[P^'fipV'\Px^>^A^^^ 


=4l 


b.fPxr^iPx'fi^'^'pf =\[Pxi,trbPx'ihb^'PZ^ + o[Pxipf^,A"bPx'fipV'\Px^>^A^^- 


(A16) 


(A17) 

(A18a) 

(A18b) 

(A18c) 

(Al9a) 

(A19b) 


Additional identities are found by permutation of the tensor product structures ri(g)r 2 g)r 3 appearing on both sides 
of the above equations. For example, applying the permutation Fi g) F 2 g) F 3 —)■ F 2 g) Fi g) r 3 to the LHS of either 


It is for this reason that we do not consider a tensor reduction technique such as the “Greek projections” used in Ref. m that commutes 
with algebraic relations valid in D-dimensions. The Greek projections provide algebraicly consistent continuations of spin-Fierz relations 
between Dirac structures of the form to D-dimensions and for example specify b\ = —80, &2 = —12 in Eq. | |A13| |. Straightforward 

generalizations of the Greek projections can relate structures of the form Fi 0 r 2 0 Fs. However, there is no straightforward extension 
of the Greek projections for Eq. | |A12| | unless (j 0 cr operators are included in the physical basis as in Ref. m- Since the nn basis of 
interest for many BSM models includes the scalar diquark operators considered here, Greek projections do not give a useful way to 
define our one-loop evanescent counterterm basis. After choosing this one-loop evanescent basis, is fully determined and we have 
no further need to establish concrete evanescent basis conventions. 
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equation leads to a new identity with the RHS modihed byl^cr^crOtT®!®^, l(g)l(g)l and a ® a ® 1 left 
unchanged, and cr 0 cr 0 cr —)■ —a ® a ® a. Identities involving general permutations of Fi 0 r 2 ® Fa are constructed 
analogously, and in particular a ® a ® a will change sign under any permutation ofFi ( 8 )F 2 ( 8 )F 3 with odd signature. 
All other Dirac structures produced by two-loop diagrams can be related to those above and our basis structures by 
algebra valid in general D. 


3. Flavor Algebra 


A convenient basis for su(2)^ tensors is given by 


^2 A 
T T , 


T T T , 


(A20) 


where the are normalized as Pauli matrices = 2. 


After applying the a 0 a spin Fierz identity of Eq. (A12), the resulting spin-singlet diquarks no longer have their 


flavor indices contracted with one of the basis structures above. Flavor (as well as color) Fierz relations are useful in 
relating the resulting structures to the original operator basis. One-loop diagrams involving flavor singlet diquarks 
require 


^2 2 _ ^^2 2 , ^/2^A\ (^2A\ 

^adJcb ~ 2 ^ ^ T jcd- 


(A21) 


Similarly, flavor vector-singlet structures require 


^2 1^2 A\ _^^2 1^2 A\ . ^ (2 B\ (2AB\ 
^adV ^ )cb — 2^abV ^ Jcd + 2 ^ )ab{T T T )cd- 


(A22) 


Finally, flavor vector-vector structures require 


l-r-‘2^A\ /^2B\ _^Sf^2^A\ (^2B\ ,(2B\ (^2A\ ,-ABC\l2n\ 2 2/2C\ 1 

(j ^ jadV T jcb ^ i ^ )ab\T T )cd + (T T )ab(T T )cd -k [(T T jabT^cd ^ab\^ ^ )cd\ 

cAB [2 2 (J 2 C\ tJ 2 C\ 


I JjAB r 2 2 1^2 C\ (^2C\ 1 1 

d” ® \jab^cd ^ )ab\T ^ )cd\ / • 


(A23) 


Eq. (A23l implies in particular that symmetric traceless tensors are Fierz self-conjugate. 


Appendix B: Two-Loop Integrals 

When evaluating simple Feynman diagrams, one can often perform Dirac “numerator algebra” that reduces the 
diagram to a simple product of a Dirac structure times a scalar integral. When evaluating diagrams with gluon 
propagators connecting quarks in separate spin-singlet diquarks, this is not possible. One is forced to work with 
tensor integrals that contain free Lorentz indices contracted with structures such as cr^j/ 0 Upr- In this case, tensor 
reduction techniques similar to those described in Appendix can be used to express tensor integrals in terms of 
linear combinations of scalar integrals. In our calculation of the diagrams of Fig. the complete set of two-loop 
tensor integrals appearing in these diagrams was organized according to the propagator powers and loop-momentum 
vectors appearing. Each tensor integral was then expressed as a linear combination of basic tensors and two-loop 
scalar integrals by tensor reduction techniques. The two-loop scalar integrals were recursively evaluated as described 
below and the results tabled for use in tensor integral evaluation. Computer algebra was essential for this process and 
performed using Mathematica scripts written by the authors. 

There exists a vast literature on evaluation of multi-loop tensor and scalar integrals. References to reviews and 
original literature are given below, and it should be emphasized that none of the techniques reviewed in this appendix 
are novel. Our aim is simply to consolidate known techniques needed for two-loop anomalous dimension calculations 
without detailing the additional complications and generalizations needed for more complex higher-order calculations. 


1. Two-Loop Scalar Integrals 

We are only concerned here with calculating the 1/e single- and double-pole pieces of two-loop diagrams. This 
allows for substantial simplifications. In particular, external momenta can be freely chosen diagram by diagram. To 
see this, note that in a renormalizable theory, these pole pieces can be at most polynomial in external momenta. 
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After factoring out possible overall dimensionful factors common to all diagrams, this means the pole pieces are 
independent of external momenta. This holds for individual diagrams as long as they contain no subdivergences, and 
therefore for general two-loop diagrams as long as one-loop counterterm diagrams cancelling all subdivergences are 
included ITT]. We may therefore freely choose a different momentum routing convenient for each two-loop diagram 
under consideration as long as the same routing is used in all corresponding one-loop counterterm diagrams. 

The only caveat to this statement is that the choice of external momentum routing must not introduce IR diver¬ 
gences. For instance, if in a massless theory one sets all external momenta to zero then all integrals vanish identically 
in dimensional regularization. This means IR divergences have been introduced that are regulated as 1/e poles and 
cancel all of the original UV divergences [5S] . We are interested in the counterterms needed to cancel UV divergences 
only, and so we must use care to choose momentum routings free of IR divergences. See Ref. [96] for a detailed review 
of this “infrared rearrangement” trick; for our purposes it is enough to note that IR divergences can be found and 
avoided through standard power counting arguments used to determine a diagram’s degree of UV divergence ITT]. 

For all diagrams in Fig. I^ a momentum routing can be chosen so that the only scalar integrals appearing are of 
the form 

/ d^kd^q 1 

(p + (p + ^)2«.^2n3g2„4(fc - q)2n, ’ 


where p is an arbitrary external momenta that serves as an IR regulator and we are suppressing omnipresent ie terms 
in factors such as -I- ie)^^. If one of the propagator factors does not appear (one rn equals zero), then the two 
loop integral can be expressed as a product of one loop integrals. The second loop includes non-integer propagator 
powers, but can still be evaluated through the textbook formula 


I{a,P) = J 


d^k 


(27r)-® {p -I- A;)2“fc2/3 


iP^) 


,2\2-a-/3-e 


T{a + l3-2 + £)r(2 - a - e)r(2 - /3 - e) 


(47r)2-'^ [ r(4 - a -/? - 2e)F(a)r(/3) 

Scalar two-loop integrals with (at least) one zero argument are given by 


r(ni,n2,n3,n4,0) = /(ns, ni)/(n4, 712 ), 

r(0, n2, ns, n4, ns) = (p^)”3+"'*“^+®/(ns, n5)J(n3 -|- n4 -|- ns - 2 -|- e, n2), 
r(ni, n 2 , 0, n4, ns) = (p^)”’’+"''“^+®/(ni, ns)/(n4, ni -I- n 2 -I- ns - 2 -|- e). 


(B2) 


(B3) 


The cases of n 2 = 0 and n 4 = 0 can be found by the (ni, ns) O (n 2 , n 4 ) symmetry of T(ni, n 2 , ns, n 4 , ns). 

This leaves the case of non-vanishing ni, • • • ,ns. This case can be evaluated recursively through the “integration 
by parts” technique of Refs. ETIIM], see Ref. IMI for a review. The starting point for this technique is the observation 
that there are no surface terms when integrating a total derivative in dimensional regularization that is 


0 = /i 


4s 


d^kd^q 

(27r)2-D 



a^{k,q,p) 


(B4) 


where a^{k, q,p) is an arbitrary vector that may depend on loop and external momenta. Useful identities are generated 
by taking to be a loop momentum vector times the integrand of Eq. (B11. Consider in particular 


d 

dqi^ 


(k - qy 


{p + fc)2"i (p -\- qj 2 ri 2 ]^ 2 n 3 q2ni _ q,^2n6 


2 n 2 {k - q) ■ {p + q) 2nyk - q) ■ q 

- ? —rW 2 - 2 -^ 2ns 

[p + qy q^ 


1 


(B5) 


{p + ky'^^ip -f g)2"2/i;2n3q2n4^;j _ g^2n5 

Next, re-write all scalar products appearing in Eq. (B51 in terms of linear combinations of p^ and denominator factors, 
for instance 


2{k -q)-{p + q) = {k+ pf - {k - qf - {p + qf, 
2{k — q) ■ q = k^ — {k — qY — q^. 


(B6) 


This allows us to express Eq. (B4) as 


0 = [2ns-I-n 2 -I-n 4 —/?-I-n22+(5 —1 )-l-n44+(5 —3 )] r(ni, n 2 , ns, n 4 , ns), 


(B7) 
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where we define 


l=‘=T(ni,n 2 ,n 3 ,n 4 ,n 5 ) = T{ni ± 1 ,712,ria,714,ris), 


(B 8 ) 


etc. This identity is sufiicient to derive a recursive solution for T(r 7 i, 712 , n. 3 , 714 , 715 ) with all rii non-zero, 

1 


T(77i,772,773,774,775) = 


H — 772 — 714 — 2775 


[t722+(5 - 1 )-|-7744+(5 - 3 )] T(77i, 772,773,774,775). (B9) 


This recursion terminates when each integral on the RHS has at least one Ui zero and the base case Eq. (B3| can 
be applied. Many other integration by parts identities and more powerful recursive algorithms can be constructed 
but are not needed for the calculation at hand. For further discussions of more general one-loop scalar integrals see 
Refs. |100L 1101] . For further discussions of two-loop scalar integral evaluation see Refs. [10211103] and the review 
Ref. IMl- 


2. Two-Loop Tensor Integrals 

Two-loop tensor integrals can be expressed in terms of scalar integrals through tensor reduction techniques. Consider 
for example the rank two integral 


Tl^{nx,n2,nz,ni,n^) = J 


4 ^ f d^kd^q 


kf tk^ 


(27r)2-® {p + k)^^^{p -I- 

By Lorentz invariance, the integral can be expressed as a linear combination 

~ '^s 9^il' + T^—p^pi,. 


Contracting both sides with these same tensors gives the system of equations 

A-2e l\ Y \ 

V 1 V 


(BIO) 


(Bll) 


(B12) 


The contractions of the RHS can be reduced to linear combinations of scalar integrals by re-writing tensor products 
in terms of differences of propagator factors as before. 


9^''Tl^{'ni,n2,nz,n4:,n^) = 3 T( 77 i, 772 , 773 , 774 , 775 ) 

1 1 r - - 21 1 (B13) 

-^P^p''T^^{ni,n2,n2„ni,n5) = ^p^^ [1 -3 - p J T^(t7i, 772 , 773 , 774 , 775 ). 

p zp 

This final formula does not apply to the cases of 777 = 0 and 773 = 0. These must be treated separately, and a 
general method can be constructed by first performing a tensor reduction of a one-loop subintegral. This problem 
is systematically considered in Ref. |104j . The following recipe is sufficient for the integrals considered in this work: 
first evaluate the one-loop integral for the loop momentum that only appears in two propagators by a one-loop tensor 
reduction. A change of variables may be useful to ensure there is only one “external momentum” scale (which may be 
a linear combination of p^ and the other loop momentum) that needs to be included in the one-loop tensor reduction. 
The second integral will then be another one-loop tensor integral involving a single scale that can be readily evaluated. 
For further discussion of tensor integral reduction techniques, see Ref. |10511106| and for a review see Ref. m- 


Appendix C: Evanescent Operators 


The MS and RI-MOM renormalization schemes are fully defined by th e ren ormalization conditions of Sec. IV and 


specification of the one-loop evanescent counterterms appearing in Sec. VB 
one-loop RI-MOM matching factors 


Two-loop MS anomalous dimensions. 


and Wilson coefficients from one-loop BSM matching all separately depend on 
the basis chosen for evanescent operator counterterms. In particular, loop-level BSM matching calculations must use 
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the same evanescent basis used in this week. Our basis includes the following evanescent operators needed as one-loop 
counterterms to Qi, 

= ('4jCPRaf_,jT2ijj){tljCPRafj,,yiT2ip){tj}CPRiT2T+ip)T^^^ - 12 Qi, (Cla) 

= {'ij)CPRiT2'il^){'iljCPRa^jT2'ip){'ij^CPR(j^^iT2T+'ip)T^^^ 

+ {'ijjCPRa^^iT2ilj){'il^CPRiT2'il^){ipCPRa^^iT2T+'ilj)T^^^ - 8Qi, 

E^ = {ipCPL(Tfj,jT2ilj){ijjCPRaf,i,iT2tp){iljCPRiT2T+^)T^^^, (Clb) 

E 2 = {ipC PLiT24’){4’CPRaf_„yiT2'il^){'ipC PRai^jT2T+'ip)T^^^ 

+ {'ll}CPLU^i^iT2lp){'ll}CPRiT2llj){ll}CPRa^uiT2T+1p)T^^^ - 4(^2, 

E^ = {'ipCPLcrfj,^iT2'tlj){'ipCPLcrfj,^iT2ip){'tljCPRiT2T+'il})T^^^ - 12 Qi, (Clc) 

E^ = {ll}CPLiT2'ip){'tpCPL(Tf,uiT2'lp){'ipCPRa^^iT2T+ll})T^^^ 

+ {lpCPL(Ti,jT2'llj){lljCPLiT2'lp){'ipCPRai,jT2T+ll;)T^^^, 

Ei = {'tl;CPRaf,i,iT2T3jp){'ijjCPRa^^iT2T3il;){'tjjCPRiT2T+'ijj)T^^^ (Old) 

- ^{ipCPRa^i,iT2TA'4!)('il}CPRa^,^iT2TA'4’){ipCPRiT2T+'ilj)T^^^ 

+ {iljCPRiT2T3ilj){'il:C PRa^^iT2T3tlj){iljCPRa^^iT2T+ip)T^^^ 

- ^{ilCPRiT2TA^){'il^CPRa^^iT2TA'li)){'ll)CPR(J^^iT2T+ll))T^^^ 

+ {'il}CPRa^^iT2T3il}){'ipCPRiT2T3'4>){%l2CPRa^^iT2T+ip)T^^^ 

- ^{'il^CPRa^,^,iT2T3ij;){ipCPRiT2T3'il}){'il)CPRa^^,iT2T+ip)T^^^ - I2Q4, 

E3 = {ll}CPRai,jT2T_'lp){'ipCPL(Tf,uiT2T+'ip){tpCPLiT2T+ll})T^^^ (Cle) 

-b {ll}CPRiT2T_1p){'4>CPL(T^yiT2T+'4>){ll^CPL(T^_i^iT2T+ll^)T^^^ 

+ { 4 lCPRa^„iT 2 T_i^){ll}CPLiT 2 T+' 4 j){'tl)CPLCrf^t,iT 2 T+'tl))T^^^ - AQ 5 , 

El = ^{ipCPRa^^iT2TAip){ipCPRa^i,iT2TA'ilj){'il}CPRiT2T+'il})T^^^ (Clf) 

+ ^{'^CPRiT2TA'll’){'lpCPRa^,^iT2TAll^){'ll^CPRa^jT2T+'lp)T^^^ 

+ ^{'iljCPRa^^iT2TA'il}){iliCPRiT2TA'ili){'il}CPRaf,^iT2T+'ilj)T^^^ + ^Qi, 

E3 = ^{'lljCPLa^„iT 2 TA'll^){i’CPL(Tfj,,yiT 2 TAlp){'tpCPRiT 2 T+-lp)T^^^ (Clg) 

+ ^{,i’CPLiT2TAi)){lpCPLCr^uiT2TA'il^){,i’CPRa^yiT2T+'ll})T^^^ 

+ ^{'^CPL<7^jT2TA'll’){'ipCPLiT2TAll’){'ll’CPRa^,^,iT2T+1p)T^^^ + AQ3. 


Above we have grouped evanescent operators that make similar contributions one-loop counterterm diagrams, see 
Tables |V)|vg 

The coefficients of Qi appearing above are determined by demanding that the RHS vanish in I? = 4. We have 
calculated them using two independent methods for verification: first by pen and paper application of the spin- 
color-flavor Fierz relations derived in Appendix]^ and second by automated Mathematica application of the operator 
projectors of Eq. 21 to explicit vertex functions constructed for each structure. It is straightforward to verify that all 
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other structures produced in amplitudes for d = 1-3 vanish by quark exchange antisymmetry. 

When constructing one-loop counterterm diagrams for two-loop diagrams d = 4-6, 16-24, and 32-45 containing d = 
1-3 subdiagrams, it is useful to employ a different evanescent operator basis Ej. The E'^ basis is defined by demanding 
that the prescription 


{8[CP^r^[CP^r^ - 4[CP^r^[CP^r^) , (C2) 

always provides valid operator identities in general D when E'^ operators are included. Applying this prescription to 
the amplitudes for d = 1-3 provides an explicit construction of the E'j operators 

Ef = {'ijjC PRafj^i,iT 2 ')p){^C PRa^i,iT 2 'ip){'iljC PRiT 2 T+')p)T^^^ (C3a) 

-8{i:^[CPRr^iT2^P){^^[CPRrf^iT2^^){r[CPR]'^‘^iT2T+i;‘^)T^^^ 

= El + 6{Qi — Qi) 

E^' = {il}CPL(Jfj,,yiT2ilj){'tpCPLafj,jT2ip){iljCPRiT2T+'if;)T^^^ (C3b) 

-8{^^[CPLr^lT2i;^){ij'^[CPLr^iT2^P^){i:^[CPRr<tT2T+^i)T^^^ 

= E- + 6{Q3 - Qs) 

E'l = ^{iljCPRai,^iT 2 TAip){tljCPRa^i,iT 2 TAilj){'il}CPRiT 2 T+'il})T^^^ (C3c) 

+ ^{dPCPRiT2TA'll^){i^CPRaf,^iT2TA'ip){'ll’CPR(Ji,^iT2T+'lp)T^^^ 

+ ^idPC PRa i,jT2TA'ip){ll’C PRiT2TA'ip){'ll’C Pr(J i,^iT2T+'lp)T^^^ 

- 5(^“ [CPRrPiT2TA^^){iP^[CPRV^iT2rAi>^){r[CP r] '"'zr2r+ 

= Al + y(Ql-Ql) 

E'a = ^{lljCPL(Tf,^iT2TA'll^)ill’CPL(7t,uiT2TA'fp){ll’CPRiT2T+ll;)T^^^ (C3d) 

- ^{^^[CPLr^T2TAi’^){'^'^[CPLr^iT2TA^^){'^^[CPR]^^lT2T+^‘^)T^^^ 

= E 3 + 2 {Q 3 — Q 3 ) 


with all other Ej equal to the corresponding Ej. The Aj basis is convenient for two-loop diagram evaluation, but is 
cumbersome for RG evolution because it includes one-loop mixing between Qi and Qi and between Q 3 and Q 3 . 

After evaluating two-loop diagrams in the E'j basis to determine the loop coefficients defined in Eq. (611, 

a change of basis to the Ej basis can be performed to recover the coefficients [Ttot]// directly appeaing in the 


anomalous dimension formula Eq. (57). Reading off the coefficients from the entries for d = 1-3 in Tables 
the necessary change of basis formula to recover ^ from are found to be 


V VI 


l(2).l ^ r^/ 


,(2)^1 _ ^,^^(0) _^r(0)>, 
'tothi 9 V 1 '1 J 

,(o) Ao)^ 


[AtotJii 

\r l(2).l _ rr/ l(2).l _ 3/„(0) _ AO)', 

[zjtotJss — l-z^totJas 2^3 ^3 ) 

\J l( 2 ).l _ rr/ l( 2 ).l _ 5/ ( 0 ) _ AO)', 

r? 1(2)4 _ rr/ 1(2)4 _ 3 /( 0 ) _ AO)', 

[rjtotJas — [J-'tothd. 2^3 ^3 )■ 


(C4) 


One-loop matching results give = 7 = . 
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Appendix D: Diagram Results 

After choosing a convenient momentum routing, each two-loop diagram in Fig. and associated counterterm 
diagrams are evaluated in terms of a tensor integral contracted with a Dirac tensor, as discussed in Appendix]^ The 
amplitude for each diagram is represented by a complicated combination of color factors and Dirac structures that 
can be simplified with the tensor reduction techniques of Appendix]^ We find it convenient to first perform a Dirac 
tensor reduction using the prescription of Eq. (29) with oi = 02 = 0 (that is, working in the E'j basis) that allows 
us to express the diagrams in a given class as individual color factors times a common combination of Dirac basis 
structures. These Dirac structures are shown in Table [Ill| We then perform color tensor reductions of the color factors 


of Table IV This allows each diagram amplitude to be expressed as a combined spin-color tensor. While there are 
25 distinct spin-color tensors that can be built from the basis tensors of Appendix]^ most vanish by quark exchange 
antisymmetry when contracted with external quark fields and flavor tensors to form operator corrections. 

When contracted with {ipiT2'ip){4’'iT2ip)i'ipiT2TA'^l^), the only spin-color tensors that give non-vanishing contributions 


are 


(l(g)l(g)l)T 


AAS 


(C 


i)r' 


sss 


(1 (g) cr (g) a)T 


ASA 


(cr (g) 1 (g) a)T 


AAS 


(Dl) 


The overall contribution from each diagram class to is represented by a linear combination of these four spin-color 
tensors in Table [v| Analogously, when contracted with (V’*T2'7gi'0)('*/'*'^2TB'i/’)('0*T2TcV^); the only spin-color tensors 
that give non-vanishing contributions are 


(1 (g) 1 (g) i)r' 


sss 


(c 


1)T 


AAS 


(1 (g) tr (g) a)T 


SAA 


(cr g) 1 (g) (j)T 


ASA 


(c 


■ (g) a)T 


AAA 


(D2) 


The overall contribution from each diagram class to Alf is represented by a linear combination of these five spin-color 
tensors in Table ED 

After adding the amplitudes Aif and Alf from each diagram class, we consider each Qi independently form 
operator corrections by Eq. ( |25[ ). The resulting operator corrections can be expressed as a simple multiple of Qj in 
D = A through application of either the operator projectors Eq. (22) or the relations of Appendix [C] and 


- (V'C PRiT2TA(T^ui^) {'li)CPR(J^piT2TA'4>) {tpC PRa^piT2T+'ll^)T 
^ {ipcPRiT2TA'ip) (iIjCPR iT2TAllj) {'lpCPRiT2T+'lp)T^^^ , 


AAA 


(D3) 


and 


i{'^CPRiT2T[3ap^4j){^CPR(7ppiT2T3'ip){tjjCPRa^piT2T+y'il;)T 


AAA 


- -i{lpCPRiT2T{A(^p^lll){'ll}CPRappiT2TA'li)){ll}CPR(J^piT2T+}ll))T 


AAA 


D=4 


= 4 


{'tpCPRiT2T[3<^p^1p){tj}CPRiT2T3'4j){ljjCPRiT2T+y4’)T 


iSSS 


(D4) 


sss 


--{'lpCPRiT2T{A'4’){'4’CPRiT2TAi’){'li^CPRiT2T+}1p)T 

The result is an operator correction proportional to Qj plus irrelevant evanescent contributions. The 1/e pole 
coefficient of this amplitude represents A [LJoJ jj^’^/e. After including the change of evanescent basis factors 


in Eq. (|C4[), is immediately given by Eq. (|57[). 
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d 

Nd 

1 (g) 1 ® 1 

(T (g) (7 (g) 1 

1 (g) (7 ® (T 

CT (g) 1 (g) (7 

(7 g) (7 (g) (7 



1/e^ 

1/e 

1/e^ 

lie 

1/e^ 

1/c 

1/e^ 

1/e 

1/e^ 

1/e 

1 

3 

- 

-4 

- 

0 

- 

0 

- 

0 

- 

0 

2 

6 

- 

-1 

- 

1/4 

- 
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- 

0 

- 

0 

- 

0 

4 

3 

-8 

8 

0 

0 

0 

0 

0 

0 

0 

0 

5 

6 

(-1 - 3dx)/2 

(5 + 3dx)/4 

1/2 

-1/2 

0 

0 

0 

0 

0 

0 

6 

6 

(-1 - 3dx)/2 

(5 + 10dx)/4 

-1/2 

15/16 

0 

0 

0 

0 

0 

0 

7 

3 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

8 

6 

0 

-2 + 3(5x 

0 

1/4 

0 

0 

0 

0 

0 

0 

9 

6 

0 

-2 + 3(5x 

0 

-1/4 

0 

0 

0 

0 

0 

0 

10 

6 

2 

2 

0 

0 

0 

0 

0 

0 

0 

0 

11 

12 

1/2 

0 

-1/8 

-5/16 

0 

0 

0 

0 

0 

0 

12 

12 

1/2 

0 

1/8 

5/16 

0 

0 

0 

0 

0 

0 

13 

6 

-2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

14 

12 

-1/2 

0 

1/8 

1/16 

0 

0 

0 

0 

0 

0 

15 

12 

-1/2 

0 

-1/8 

-1/16 

0 

0 

0 

0 

0 

0 

16 

12 

-2 

1 + 7SJA 

0 

-1/16 

0 

0 

0 

0 

0 

0 

17 

12 

-2 

-1 

1/2 

1/2 

0 

0 

0 

0 

0 

0 

18 

12 

-2 

1 - 7dx/4 

0 

1/16 

0 

0 

0 

0 

0 

0 

19 

12 

-2 

-1 

-1/2 

-1/2 

0 

0 

0 

0 

0 

0 

20 

12 

-1/2 + 3(5x/2 

(3 - 3dx)/4 

-1/4 

1/8 

0 

0 

0 

0 

0 

0 

21 

12 

-1/2 + 3<5x/2 

(3 - 10dx)/4 

1/4 

-9/16 

0 

0 

0 

0 

0 

0 

22 

3 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

23 

3 

-1 - 3dx 

0 

1 

0 

0 

0 

0 

0 

0 

0 

24 

3 

-1 - 3dx 

7SJ2 

-1 

7/8 

0 

0 

0 

0 

0 

0 

25 
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-6 

5 

0 

0 

0 

0 

0 

0 

0 

0 

26 

12 

-3/2 

1/2 

3/8 

1/16 

0 

0 

0 

0 

0 

0 

27 

12 

-3/2 

1/2 

-3/8 

-1/16 

0 

0 

0 

0 

0 

0 

28 

12 

0 

0 

0 

3/4 

0 

0 

0 

0 

0 

0 

29 

3 

15/2 - Nf 

-13 + 4Ai//3 

0 

0 

0 

0 

0 

0 

0 

0 

30 

6 

0 

0 

-5/8 +Ar//12 

17/48-Ary/72 

0 

0 

0 

0 

0 

0 

31 
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0 

0 

5/8 - Nf/12 

-17/48 +Ary/72 

0 

0 

0 

0 

0 

0 

32 

6 

-4 

0 

1 

0 

0 

0 

0 

0 

0 

0 

33 

6 

-4 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

34 

6 

-1/2 

0 

1/8 

1/16 

-Ax/8 

-1/16 +Ax/16 

1/8 

-1/16 

-1/8 

0 

35 

6 

-1/2 

0 

1/8 

1/16 

Ax/8 

1/16-A^/ie 

-1/8 

1/16 

1/8 

0 

36 

6 

-1/2 

0 

1/8 

1/16 

-Ax/8 

-1/16 +Ax/16 

1/8 

-1/16 

1/8 

0 

37 

6 

-1/2 

0 

1/8 

1/16 

Ax/8 

1/16-Ax/16 

-1/8 

1/16 

-1/8 

0 

38 

6 

-1/2 

0 

-1/8 

-1/16 

-Ax/8 

-1/16 +Ax/16 

-1/8 

1/16 

-1/8 

0 

39 

6 

-1/2 

0 

-1/8 

-1/16 

Ax/8 

1/16-Ax/16 

1/8 

-1/16 

1/8 

0 

40 

6 

-1/2 

0 

-1/8 

-1/16 

-Ax/8 

-1/16 +Ax/16 

-1/8 

1/16 

1/8 

0 

41 

6 

-1/2 

0 

-1/8 

-1/16 

Ax/8 

1/16-Ax/16 

1/8 

-1/16 

-1/8 

0 

42 

6 

-1 

0 

1/4 

0 

-Ax/4 

0 

1/4 

0 

-1/4 

0 

43 

6 

-1 

0 

1/4 

0 

Ax/4 

0 

-1/4 

0 

1/4 

0 

44 

6 

-1 

0 

-1/4 

0 

Ax/4 

0 

1/4 

0 

1/4 

0 

45 

6 

-1 

0 

-1/4 

0 

-Ax/4 

0 

-1/4 

0 

-1/4 

0 

46 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3/8 


TABLE III: NDR 1/e pole structure of the diagram amplitudes in Feynman gauge without color factors, d labels the diagrams 
classes of Fig. Nd is the number of diagrams within class d, 5^ = dxiXa = ^XiX 2 ^X 2 X^ ■ Evanescent counterterms in 

the Ej basis defined by applying Eq. (|60[) to diverent subdiagrams are included in these results. 
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d 

Color Factor 

1 

1 ■ 
2 . 

rp ^ T' 

4 (ji)(fcZ){mn} 3 4 (ij)(fcZ){mn} 


2 

1 • 
2 . 

rp ^ T' 

4 (iZ)(fej){mn} 3 4 (ij)(fcZ){mn} 


3 

1 ■ 

2 . 

T(ik)(jl){mn} — 3 T(ij)(kl){mn} 


4 

1 1 2 rri 1 

4 L 3 -‘-Ui)(kl){mn} + g J-(ij)(kl){mn}\ 

5 

1 \ 2 rri I 10 ^ 1 

4 L + g J-(ij)(kl){mn}\ 

6 

1 1 2 n-\ 1 

4 L 3 -‘-(ik)(jl){mn} + g J-(ij)(kl){mn}\ 

7 

4 + 9'^(ij)(kl){mn}] 

8 

4 [^'^(iO(kj){mn} + gT(ij) (kl) (mn}] 

9 

4 [^'^(ik)Ul){mn} + g Ffij)(fei) 1™"}] 

10 

1 

12 

rp ^ 'T' 

4 (ji)(fcZ){mn} 3 4 (ij)(fcZ){mn 


11 

1 

12 

rp ^ 'T' 

4 (iZ)(Ze}){mn} 3 4 (ij)(fcZ){mn 


12 

1 

12 

T(ik){jl){7nn} ~ '^T(^ij){kl){mn' 


13 

3 ['^Ui)(kl){mn.} “ 3 Fcij)(fci){mn}] 

14 

3 ['^(il)(kj){mn.} “ 3 F(ig)(*;;){m„}] 

15 

3 ['^(ik)Ol){mn.} “ 3 F(ig) (*;() {^n}] 

16 

1 

4 

T(jl)(ki){mn} 3 T(ji)(kl){mn} 3 T(il)(kj){mn} + g F(ig)(fcj) 


17 

1 

4 

<T-T 1 rj-\ 1 rp 1 1 rji 

{li){kj){mn} 3 (ji)(fcZ){mn} 3 (iZ)(fcj){mn} "r g(ij)(fcZ){mn} 


18 

1 

4 

'rp 1 rp 1 ry-T . 1 rp 

{jk){il){mn} 3 (ji)(fcZ){mn} 3 (ifc)(jZ){mn} "r g(ij)(fcZ){mn} 


19 

1 

4 

<T-T 1 rp 1 ry-T . 1 rp 

{ki){jl){mn} 339 


20 

1 

4 

J-{ik){lj){mn} 3 (zZc)(jZ){mn} 3 (iZ)(fcj){mn} "r g 4 (j y)(fc^){Tnn} 


21 

1 

4 

T(il)(jk){m.n.} — 3 T(ik)(jl){mn.} ~ 3 T(i;) (fcg) + g r(ig) (fcl) {mn} 


22 

1 

4 

T(ji)(lk){mn} ~ 3 T(ij)(lk){mn} ~ 3 T(ji)(kl){mn} + g F(ig)(fcj) 


23 

1 

4 

T(kl)(ij){mn) ~ 3 T(il)(kj){mn} ~ 3 T(kj)(il){mn} + g F(ig)(fcj) 


24 

1 

4 

T(lk)(ji){mn) ~ 3 T(lj)(ki){m.n} ~ 3 T(ik) (j 1) {mn} + g F(ig) (fcj) 


25 

3 

4 

T(^ji)(kl){mn} 3 T(ij)(kl){mn} 


26 

3 

4 

T(il)(kj){mn} ~ ~^T(ij)(kl){mn} 


27 

3 

4 

rp ^ T' 

(ik)(jl){mn} (ij)(kl){mn} 


28 

1 

4 

'^{jl)(ki){mn} ^(Zi) (fcj) {mn} ] 

29 

1 ■ 

2 . 

rp ^ T' 

4 (ji)(fcZ){mn} 3 4 (ij)(A;Z){mn} 


30 

1 • 
2 . 

rp ^ T' 

4 (iZ)(fcj){mn} 3 4 (ij)(fcZ){mn} 


31 

1 ■ 

2 . 

T(ik)(jl){mn.} ~ 3 T(ij)(kl){mn.} 


32 

1 

4 

T(il)(kj){nm} ~ 3 T(il)(kj){mn} ~ 3 T(ij)(kl){nm} + g F(ig)(fcj) 


33 

1 

4 

T(lj)(ki){nm} ~ 3 T(lj)(ki){mn} ~ 3 T(ij)(kl){nm} + g F(ig)(fcj) 


34 

1 

4 

T(^in)(kj){ml} ~ 3 F(ij) (fcj) {rnn} “ 3 F(i„) (fcj) + g (j.;) 


35 

1 

4 

/yn 1 /y-i 1 /y-i 1 rp 

J- (im){kj){ln} ~ J-(il)(kj){mn} ~ J-(im){kl){jn} ~r g J-(ij){kl){mn} 


36 

1 

4 

rp 1 rp 1 rp . 1 rp 

{mj){il){kn} 3 (Zcj) (iZ){mn} 3 4 (mj) (fcZ){in} "r g 4 (ij )(fcZ){mn} 


37 

1 

4 

rp 1 rp 1 ry-T . 1 rp 

4 (nj)(iZ){mfc} 3 (Zcj)(*0 3 (^i)(^0 9 (^0 


38 

1 

4 

rp 1 rp 1 rp . 1 rp 

4 (im)(}Z){fcn} 3 4 (iZc)(jZ){mn} 3 4 (irn)(fcZ){}n} "r g 4 (j y)(fc^){Tnn} 


39 

1 

4 

ry-T 1 ry-i 1 rp 1 rp 

4 ( 2 n)(jZ){mfc} 3 4 (iZc) (jZ){mn} 3 4 (in) (fcZ){m}} “1” g 4 (ij) (fcZ){mn} 


40 

1 

4 

T(nj)(ki){ml) ~ 3 T(lj)(ki){mn} ~ 3 T(nj)(kl){mi} + g T(ig) (fcj) 


41 

1 

4 

F(mi){*;i){ln} “ 3T(lj)(ki){mn} ~ 3T(mj)(kl){in} + g F(ig)(fcj) 


42 

1 

4 

T(ml)(kj){in} ~ 3 F(ij) (fcg) {rnn} “ 3T(mj){kl){in} + g F(ig) (fcj) 


43 

1 

4 

T(^nl)(kj){mi} ~ 3 F(ij) (fcg) {rnn} “ 3T(„j)(kl){mi} + g F(ig) (fcj) 


44 

1 

4 

rp 1 /y-i 1 /y-i 1 ry-i 

4 (mfc)(jZ){in} ~ 3 4 (ifc)(jZ){mn} ~ 3 4 (mj)(fcZ){in} "T g 4 (jj ) (fcZ){mn} 


45 

1 

4 

ry-T 1 rp 1 rp . 1 rp 

J-{nk){jl){mi} 3 4 (ifc)(jZ){mn} 3 4 (nj)(fcZ){mi} "r g 4 (ij ) (fcZ){mn} 


46 

4 [^(®")( fef ){™0 


TABLE IV: Single diagram co lor f actors corresponding to the explicit diagrams in Fig.[^ This table alleviates potential sign and 


coefficient ambiguity in Table III due to choice of color terms factored out. T^ij)(ki){mn} represents T^j][ki]{mn} o'" 


depending on which operator is inserted in the diagram. 
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d 

(1 (g) 1 (g) 

(o- (g) (T (g) 

(1 (g) (7 (g) + (ct (g) 1 (g) 


l/e^ 

1/e 

l/e^ 

1/e 

l/e= 

1/e 

1 

- 

4 

- 

0 

- 

0 

2 

- 

3/2 

- 

-1/8 

- 

-1/8 

3 

- 

3/2 

- 

-1/8 

- 

-1/8 

4 

-8 

8 

0 

0 

0 

0 

5 

(-11 - - 135^)/12 

(55 -b 75a + 135a)/24 

1/12 

-1/12 

1/12 

-1/12 

6 

(-11 - - 135^)/12 

(165 -b 705a + 1305a)/72 

1/12 

-5/32 

1/12 

-5/32 

7 

0 

1 

0 

0 

0 

0 

8 

0 

(10-b235i - 105\)/12 

0 

-7/48 

0 

-7/48 

9 

0 

(10-b 235 a - 195a)/12 

0 

-7/48 

0 

-7/48 

10 

2/3 

2/3 

0 

0 

0 

0 

11 

1/4 

0 

-1/48 

-5/96 

-1/48 

-5/96 

12 

1/4 

0 

-1/48 

-5/96 

-1/48 

-5/96 

13 

16/3 

-8/3 

0 

0 

0 

0 

14 

2 

0 

-1/6 

-1/12 

-1/6 

-1/12 

15 

2 

0 

-1/6 

-1/12 

-1/6 

-1/12 

16 

-2/3 

(24 - 285i -b 355i)/72 

0 

1/48 

0 

-1/96 

17 

-2/3 

-1/3 

1/3 

1/3 

1/12 

1/12 

18 

-2/3 

(24-b 285 a - 355i)/72 

0 

1/48 

0 

-1/96 

19 

-2/3 

-1/3 

1/3 

1/3 

1/12 

1/12 

20 

(-1 - 13<5\-h 8<5i)/12 

(3 -b 135i - 85i)/24 

1/8 

-1/16 

0 

0 

21 

(-1 - 13(5ji -b 85i)/12 

(9 -b 1305i - m \)/72 

1/8 

-9/32 

0 

0 

22 

0 

0 

0 

0 

0 

0 

23 

(-5 - 75 \ - 45i) 

0 

1/12 

0 

1/3 

0 

24 

(-5 - 75 \ - 45i) 

(495i -b 285i)/72 

1/12 

-7/96 

1/3 

-7/24 

25 

-18 

15 

0 

0 

0 

0 

26 

-27/4 

9/4 

9/16 

3/32 

9/16 

3/32 

27 

-27/4 

9/4 

9/16 

3/32 

9/16 

3/32 

28 

0 

0 

0 

-3/4 

0 

0 

29 

-15/2-b TV/ 

CO 

1 

CO 

0 

0 

0 

0 

30 

0 

0 

5/16-iV//24 

-17/96+ iV//144 

5/16- Ar//24 

-17/96 +A//144 

31 

0 

0 

5/16-iV//24 

-17/96+ iV//144 

5/16- 77^/24 

-17/96 +Ar//144 

32 

-14/3 

0 

-1/6 

0 

1/3 

0 

33 

-14/3 

0 

-1/6 

0 

1/3 

0 

34 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

35 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

36 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

37 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

38 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

39 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

40 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

41 

1/48 

0 

(10- 3Ax)/192 

(-1 + Ax)/128 

(10-3Ax)/192 

(-1 + Ax)/128 

42 

-11/24 

0 

(10-3Ax)/96 

0 

(-2-9Ax)/96 

0 

43 

-11/24 

0 

(10-3Ax)/96 

0 

(-2-9Ax)/96 

0 

44 

-11/24 

0 

(10-3Ax)/96 

0 

(-2-9Ax)/96 

0 

45 

-11/24 

0 

(10-3Ax)/96 

0 

(-2-9Ax)/96 

0 

46 

0 

0 

0 

0 

0 

0 


TABLE V: Pole structure for operators Qi, Q 2 , Qs in the Ej evanescent basis. As before d labels the diagram class. 

Unlike Table [Till the remaining columns include the total contribution from all diagrams in the class. Only spin-color tensors 
with index symmetries appropriate for these operators are shown. 5 \ = 5 xix 2 > = '^X2X3 + ^xixs- 
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d 

(1 (g) 1 (g) 

(a (g) a (g) 

+ (1 (g) cr (g) + (ct (g) 1 (g) 

(a (g) cr (g) 



1/e 

1/e^ 

1/e 

1/e^ 

1/e 

1 

- 

-4 

- 

0 

- 

0 

2 

- 

5/2 

- 

-3/8 

- 

0 

3 

- 

5/2 

- 

-3/8 

- 

0 

4 

- 8/3 

8/3 

0 

0 

0 

0 

5 

(-13- 135s)/12 

(65 + 135s)/24 

1/4 

-1/4 

0 

0 

6 

(-13- 135s)/12 

(195 + 1305s)/72 

1/4 

-15/32 

0 

0 

7 

0 

-11/3 

0 

0 

0 

0 

8 

0 

(38 - 195s)/12 

0 

-7/16 

0 

0 

9 

0 

(38 - 195s)/12 

0 

-7/16 

0 

0 

10 

-2/3 

-2/3 

0 

0 

0 

0 

11 

5/12 

0 

-1/16 

-5/32 

0 

0 

12 

^jvi 

0 

-1/16 

-5/32 

0 

0 

13 

-16/3 

8/3 

0 

0 

0 

0 

14 

10/3 

0 

-1/2 

-1/4 

0 

0 

15 

10/3 

0 

-1/2 

-1/4 

0 

0 

16 

10/3 

(-60 - 35(5s)/36 

0 

-1/8 

0 

0 

17 

10/3 

5/3 

-1/2 

-1/2 

0 

0 

18 

10/3 

(-60 + 35(5s)/36 

0 

-1/8 

0 

0 

19 

10/3 

5/3 

-1/2 

-1/2 

0 

0 

20 

(l-fc)/12 

(-3 + 5s)/24 

-3/8 

3/16 

0 

0 

21 

(l-fc)/12 

(-9 + 10<5s)/72 

-3/8 

27/32 

0 

0 

22 

-16/3 

0 

0 

0 

0 

0 

23 

(-13- 135s)/12 

0 

1/4 

0 

0 

0 

24 

(-13- 135s)/12 

91(5s/72 

1/4 

-7/32 

0 

0 

25 

18 

-15 

0 

0 

0 

0 

26 

-45/4 

15/4 

27/16 

9/32 

0 

0 

27 

-45/4 

15/4 

27/16 

9/32 

0 

0 

28 

0 

0 

0 

9/4 

0 

0 

29 

15/2 - Ns 

-13 + 4A‘//3 

0 

0 

0 

0 

30 

0 

0 

15/16-Ar//8 

-17/32 +A//48 

0 

0 

31 

0 

0 

15/16-iV//8 

-17/32 +A//48 

0 

0 

32 

10/3 

0 

-1/2 

0 

0 

0 

33 

10/3 

0 

-1/2 

0 

0 

0 

34 

-25/48 

0 

(10-3A;^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

35 

-25/48 

0 

(10-3A;^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

36 

-25/48 

0 

(10 - 3A^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

37 

-25/48 

0 

(10 - 3A^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

38 

-25/48 

0 

(10 - 3A^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

39 

-25/48 

0 

(10-3A^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

40 

-25/48 

0 

(10-3A^)/64 

(-3 + 3A^)/128 

-9/32 

0 

41 

-25/48 

0 

(10-3A;^)/64 

(-3 + 3A;^)/128 

-9/32 

0 

42 

-25/24 

0 

(10-3A;^)/32 

0 

-9/16 

0 

43 

-25/24 

0 

(10-3A;^)/32 

0 

-9/16 

0 

44 

-25/24 

0 

(10-3A;^)/32 

0 

-9/16 

0 

45 

-25/24 

0 

(10-3A^)/32 

0 

-9/16 

0 

46 

0 

0 

0 

0 

0 

9/4 


TABLE VI: Pole structure for T®®® operators Q 4 , Qs, Qi, Q 3 analogous to Table [v| Ss = 5xiX2 + '5x2X3 + 





























































